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^ ■ ABSTRACT 

Q^. We present a wide- field IFS survey on the nearby face-on Sbc galaxy NGC 628, 

comprising 11094 individual spectra, covering a nearly circular ficld-of-view of ~6 
arcniin in diameter, with a sampling of ~ 2.7 arcsec per spectrum in the optical wave- 
length range (3700-7000 A). This galaxy is part of the PPAK IFS Nearby Galaxies 

^ I Survey, (PINGS, Rosalcs-Ortcga et al. 2009). To our knowledge, this is the widest 

spectroscopic survey ever made in a single nearby galaxy. A detailed flux calibration 
was applied, granting a spectrophotometric accuracy of ^ 0.2 mag. The spectroscopic 
data was analysed both as a single integrated spectrum that characterises the global 
properties of the galaxy, and using each individual spectrum to determine the spa- 
tial variation of the stellar and ionized gas components. The spatial distribution of 
Q^ , the luminosity- weighted ages and metallicities of the stellar populations was analysed. 

t;;j- ■ Using typical strong emission line ratios we derived the integrated and 2D spatial 

distribution of the ionized gas, the dust content, SFR, and oxygen abundance. 

The age of the stellar populations shows a negative gradient from the inner (older) 
^^ ' to the outer (younger) regions. We found an inversion of this gradient in the central 

~ 1 kpc region, where a somewhat younger stellar population is present within a ring 
at this radius. This structure is associated with a circumnuclear star-forming region 
at ~ 500 pc, also found in similar spiral galaxies. From the study of the integrated 
and spatially resolved ionized gas we found a moderate SFR of ~2.4 M0 yr^^. The 

j^j ■ oxygen abundance shows a a clear gradient of higher metallicity values from the inner 

C^ I part to the outer part of the galaxy, with a mean value of 12 -I- log(0/II) ~ 8.7. At 

some specific regions of the galaxy, the spatially resolved distribution of the physical 
properties show some level of structure, suggesting real point-to-point variations within 
an individual H II region. Our results are consistent with an inside-out growth scheme, 
with stronger star formation at the outer regions, and with evolved stellar populations 
in the inner ones. 

Key words: techniques: spectroscopic - galaxies: individual: NGC 628 (M74) - 
galaxies: abundances - stars: formation - galaxies: ISM - galaxies: stellar content 



1 INTRODUCTION 

Galaxies in the local universe (~ 10 Mpc) are the funda- 
mental anchor points of any study of the evolution of these 

, „ , , . „ , , ^ » , . ,,. objects in cosmoloeical time scales. Therefore, it is impor- 

Based on observations collected at the Centro Astronomico His- 1 T 1 

pano Aleman (CAHA) at Calar Alto, operated jointly by the ^^""^ *» understand their mam properties, including their 

Max-Planck Institut ffir Astronomic and the Instituto de As- morphology, ionized and neutral gas content, stellar popu- 

trofisica dc Andalucia (CSIC). lations and metallicities. Due to their apparent scale- length 

t E-mail: sanchez@cefca.es they represent the perfect laboratories to study the depen- 



o 

C/3 



(N 
> 



o 



X 



S. F. Sanchez et al. 



dence of the stellar population, the star formation history, 
and star formation rate on the morphology and morphologi- 
cal substructures, the metallicity enrichment and the mech- 
anisms of metal transfer, as well as the nature of the gas 
ionization. 

Powerful constraints on theories of galactic chemical 
evolution, on the star formation history, and on the nu- 
cleosynthesis in galaxies can be derived from the accurate 
determination of the nature of the ionization, the star- 
formation rate and the chemical abundances at different lo- 
cations within a galaxy. Fortunately, these sort of studies 
can be addressed using ground-based astronomy by observ- 
ing galaxies with large apparent scale-lengths. Nearby galax- 
ies may be more easily separated into a number of different 
morphological components and several types of stellar pop- 
ulations. Given the spatial variation in the star-formation 
histories (including violent episodes in some cases), and the 
amount of dust attenuation within a galaxy, and their rela- 
tively low surface brightness, it is difficult to study the stel- 
lar populations of nearby galaxies using just their integrated 
light. More information can be extracted by studying their 
resolved properties, although it is a complex task to tackle. 

Nearby galaxies have been observed for many decades 
using many different techniques, such as multi-band optical 
and near-infrared broad-band imaging (to derive the prop- 
erties of their dominant stellar populations and population 
gradients), narrow-band imaging and multi-object observa- 
tions to derive their gas content and gas kinematics (e.g. 
Kennicutt & Hodge 1980; Belley & Roy 1992; Scowen et al. 
1996), and slit spectroscopy of the brightest HII regions 
within the galaxy (e.g. McCall et al. 1985; Zaritsky et al. 
1994; van Zee et al. 1998). More recent studies, using space- 
based observatories, have revealed new features in these ap- 
parently well known objects, such as star-forming regions at 
very large radii (e.g. GALEX Gil de Paz et al. 2007a, b) and 
obscured star-forming regions (e.g. Kennicutt et al. 2003; 
Prescott et al. 2007, SPITZER). 

Despite all these efforts, we still lack of a complete pic- 
ture of the main properties of these galaxies, especially those 
ones that can only be revealed by spectroscopic studies (like 
the nature of the ionization and/or the metal content of 
the gas). This is because previous spectroscopic studies only 
sample a very few discrete regions in these complex tar- 
gets (e.g. Roy & Walsh 1988; Kennicutt & Garnett 1996), 
and in many cases they were sampling very particular types 
of regions (HII regions, McCall et al. 1985; van Zee et al. 
1998; Castellanos et al. 2002b, a). Integrated spectra over 
large apertures are required to derive these properties in a 
more complete way, but are difficult to obtain using classical 
slit spectroscopy (although recent efforts have obtained inte- 
grated spectra of some local galaxies by adopting a drift-scan 
procedure (Moustakas & Kennicutt 2006b). Even in these 
cases, only a single integrated spectrum is derived, and the 
spatial information is lost. 

Recent studies have derived the integrated and spatially 
resolved properties of certain portions of nearby galaxies, 
by using Integral Field Spectroscopy (IPS) techniques (e.g. 
SAURON, Bacon et al. 2001; de Zeeuw et al. 2002). This 
technique allows one to obtain spatially resolved, continu- 
ously sampled spectroscopy over the field-of-view (FOV) of 
the instrument. These studies, though leading to extremely 
important results, were focused on the study of the cen- 



tral regions of galaxies, with a FOV of radius ~ 20 arc- 
sec, corresponding on average to the inner ~ 1 kpc. More 
recently, Blanc et al. (2009) observed the central region of 
M51 (~ 1.7 arcmin^) using the VIRUS-P instrument. How- 
ever, to date there has been no systematic study of a sample 
of local universe galaxies using IFS, covering a substantial 
fraction of their optical sizes. 

In order to fill this gap, we began an IFS sur- 
vey of 17 nearby (< 100 Mpc) galaxies, called PINGS 
(PPAK Integral-field-spectroscopy Nearby Galaxies Survey, 
Rosales-Ortega et al. 2010, hereafter Paper I). By their na- 
ture, most of the objects in the survey cannot be covered in a 
single pointing with IFS instruments, and a new observing- 
reduction technique had to be developed to perform accu- 
rate mosaicking of the targets. In this article we present the 
first scientific results based on one of our targets, NGC 628 
(or M74). This is a local universe galaxy at z ~ 0.00219 
(~9 Mpc), and the largest in projected angular size of the 
PINGS sample. 

In this first article of the series devoted to the IFS study 
of NGC 628, we present a study of the small and interme- 
diate scale variation in the line emission and stellar con- 
tinuum of NGC 628. We derive these properties from both 
the integrated spectrum of the galaxy and the spatially re- 
solved spectra (by means of pixel-resolved maps across the 
disk of the galaxy), and compare the results. Additionally, 
we include a description of the data acquisition and reduc- 
tion techniques, in particular on those details that are dif- 
ferent than the standard reduction of IFS data and that 
were not addressed in Paper I. The structure of the article 
is as follows: In Section 2 we give an overview of the general 
properties of NGC 628. In Section 3 we present the obser- 
vational details, including the instrument, telescope and the 
observing technique. In Section 4, we explain the reduction 
technique, describing the software packages used, and the 
sanity checks performed. In Section 5, we describe the anal- 
ysis performed on each spectrum included in our dataset 
and consider also the integrated spectrum of the galaxy. 
The spatial distributions of the different derived properties 
are also presented here. In Section 6 we discuss the results 
and summarise the main conclusions of this study. Finally, 
in Appendix A, we describe in detail the technique used to 
analyse the stellar populations in the galaxy, including sim- 
ulations to characterize its reliability. 



2 GENERAL PROPERTIES OF NGC 628 

Morphologically classified as an Sbc (Holmberg 1975), 
NGC 628 it is nearly face-on spiral, showing a typical grand- 
design structure. This galaxy has been well observed at 
a variety of wavelengths and there is abundant multi- 
wavelength ancillary data from photographic plates and 
CCD imaging in the optical (e.g. Holmberg 1975; Boroson 
1981; Shostak & van der Kruit 1984; Natah et al. 1992; 
Hoopes et al. 2001; Kennicutt et al. 2008) to SPITZER in 
the NIR (SINGS, Kennicutt et al. 2003), and GALEX in 
the UV (Gil de Paz et al. 2007a). NGC 628 is a good exam- 
ple of an isolated galaxy, previous studies have established 
that NGC 628 has not had an encounter with satellites or 
other galaxies in the last 10^ yr. (Kamphuis & Briggs 1992). 
However, it is possible that the galaxy suffered some kind 
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of interaction 1 Gyr ago, which could explain the presence 
of a large-scale oval structure (weak bar) discovered in the 
NIR (James & Seigar 1999; Seigar 2002), and the disturbed 
morphology in the north part of the galaxy. 

By means of observations in neutral hydrogen, an ellip- 
tical ring-like structure was discovered well beyond the opti- 
cal disc, at around 12 arcmin from the nucleus of the galaxy 
(Roberts 1962; Briggs et al. 1980), lying on a plane with 
~ 15° inclination with respect to the inner disc. The pres- 
ence of this warped velocity field is a puzzle, given the appar- 
ent isolation of the galaxy which would rule out its origin by 
tidal disruption. However, this feature is most likely the re- 
sult of the interaction with two large high- velocity clouds ac- 
creting onto the outer parts of the disc (Kamphuis & Briggs 
1992; Lopez-Corredoira et al. 2002; Beckman et al. 2003). 

From the morphological and dynamical point of view, 
NGC 628 displays one prominent spiral arm to the south, 
and one or several disturbed spiral arms to the north 
(although UV observations have shown spiral arms with 
a more symmetrical appearance than in the optical), 
an inner rapidly rotating disc-like structure (Daigle et al. 
2006), a CO-discovered circumnuclear ring of star forma- 
tion at ~ 2 kpc from the centre (Wakker & Adler 1995; 
James & Seigar 1999) thought to be the result of a barred 
potential (Seigar 2002), and a nuclear (nested) bar on a 
^ 100 pc scale (Laine et al. 2002). All these ingredients 
seem to suggest that the evolution of the structure in 
NGC 628 has been driven by secular evolution of the disc 
(Kormendy & Kennicutt 2004). 

On the other hand, Cornett et al. (1994) found that the 
star formation history in NGC 628 varies with galactocentric 
distance (with a young stellar population in the circumnu- 
clear region), while Natali et al. (1992) suggested a scenario 
with an inner and an outer disc with different stellar popula- 
tions, with a transition region located at 8-10 kpc from the 
centre, the same radius at which Cepa & Beckman (1990) 
found that the star formation efficiency is at its minimum, 
which they interpreted as the corotation radius for the spiral 
arms. 

More recently, Fathi et al. 2007 (hereafter Fathi07) pre- 
sented a detailed kinematic analysis of the galaxy based on 
wide-field 2D Fabry-Perot maps of Ha. They found that the 
velocity dispersion of the gas and for individual H II regions 
is practically constant at all galactocentric distances cov- 
ered by their study. This result, together with the fact that 
they were not able to distinguish true diffuse gas emission 
given that the star formation is widely distributed within the 
disc, let Fathi07 to suggest that the emission in Ha from 
the HII regions dominates any emission from the diffuse 
component. Furthermore, they confirmed the presence of a 
disc-like central structure, which they argue could have been 
built up by inflow from large galactocentric distances, sug- 
gesting that NGC 628 is in the process of forming a secular 
pseudo-bulge. 

NGC 628 has been also a classical target for optical 
spectroscopy, there is a good spectroscopic coverage of many 
of its structures, from the core (Moustakas & Kennicutt 
2006b) to many HII regions located within and beyond 
the optical disc of the galaxy (e.g. McCall et al. 1985; 
van Zee et al. 1998; Ferguson et al. 1998; Castellanos et al. 
2002b). Integral field spectroscopy of the central core (--^ 33 
arcsecx41 arcsec) of NGC 628 has been recently obtained 



using the SAURON Integral Field Unit (Cauda et al. 2006, 
2007). These studies showed that the stellar component ro- 
tates in the same direction as the gas, that the stellar ve- 
locity dispersion drops in the central zones (suggesting a 
dynamical cold inner disc), and that the H/3 distribution is 
more extended than the [O III] distribution, both suggest- 
ing a ring-like structure, confirming the 15-20 arcsec nuclear 
ring previously reported by Wakker & Adler (1995). 

In terms of studies focusing on the gas-phase of the 
galaxy, and on particular on the determination of the chem- 
ical abundance of the galaxy, previous long-slit spectroscopic 
studies have derived the abundance gradient of NGC 628 
up to relatively large galactocentric radii (~ 2R25), us- 
ing mainly empirical metallicity indicators based on the 
ratios of strong emission lines (e.g. McCall et al. 1985; 
Zaritsky et al. 1994; van Zee et al. 1998; Ferguson et al. 
1998; Castellanos et al. 2002b). These studies have shown 
a higher metallicity content in the inner part of the galaxy, 
that the slope of the gradient is constant across the range 
of galactocentric distances sampled by the different studies, 
and that the oxygen abundance decrease is relatively small. 
In the dynamical scenario previously described, this allows 
a moderate mixing of the disc material driven by the large- 
scale oval distortion (Zaritsky et al. 1994). However, these 
results have been drawn from relatively few spectroscopi- 
cally observed HII regions. 

An alternative method was developed by Belley & Roy 
(1992), who derived reddenings, H/? equivalent widths, di- 
agnostic line ratios and metallicities for 130 HII regions by 
the implementation of imaging spectrophotometry, i.e. using 
narrow-band interference filters with the bandpass centered 
on several key nebular lines. They found no trends in the 
reddening nor for the H/3 EW's with galactocentric distance. 
However they found that the excitation, and some diagnos- 
tic line ratios are strongly correlated with galactocentric ra- 
dius. They also derived an oxygen abundance gradient of 
NGC 628 based on the [OlIl]/H/3 ratio. 

In summary, NGC 628 represents an interesting object 
for a full 2D spectroscopic study, considering that it is 
the prototype of a normal, isolated, grand-design, face-on, 
nearby spiral, with an interesting morphology and partic- 
ular dynamical features. Furthermore, despite all the pre- 
vious studies towards a full understanding of NGC 628, we 
still lack a detailed knowledge of the gas chemistry, stellar 
populations and the global properties across the surface of 
this galaxy, which can be only obtained by spectroscopic 
means. The IFS data presented in this and future papers 
aim to fill this gap, by deriving and comparing the global 
and spatially resolved spectroscopic properties of NGC 628 
with previous results obtained by multi- wavelength observa- 
tions. The information provided by these 2D spectral maps 
will allow us to test, confirm, and extend the previous body 
of results from small-sample studies, and would provide a 
new path for the analysis of the two-dimensional metallicity 
structure of disks and the intrinsic dispersion in metallicity 
on this and other late-type normal spirals. 



3 OBSERVATIONS 

Observations were carried out at the 3.5m telescope of the 
Calar Alto observatory with the Potsdam Multi Aperture 
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Table 1. Log of the observations. 



Date 



Pointings 



Figure 1. B-band Digital Sky Survey image of NGC628. The 
mosaic of the PPAK pointings is shown as overlaid hexagons in- 
dicating the field-of-view of the central fibre-bundle. The identifi- 
cation number of each pointing is indicated, as listed in Table 1. 
The circle indicates the 1)25 radius (in the _B-band). The image 
is 10' X 10' and it is displayed in top-north, left-oast standard con- 
figuration. 



Spectrograph, PMAS, (Roth et al. 2005) in the PPAK mode 
(Verheijen et al. 2004; Kelz et al. 2006). The PPAK fibre 
bundle consists of 382 fibres of 2.7 arcsec diameter each (see 
Fig. 5 in Kelz et al. 2006). Of these 382 fibres, 331 (the sci- 
ence fibres) are concentrated in a single hexagonal bundle 
covering a field-of-view of 72" x 64" , with a filling factor of 
~ 65%. The sky background is sampled by 36 additional fi- 
bres, distributed in 6 bundles of 6 fibres each, distributed 
along a circle ~ 90 arcsec from the center of the instrument 
FOV. The sky-fibres are distributed among the science fi- 
bres within the pseudo-slit in order to have a good charac- 
terization of the sky; the remaining 15 fibres are used for 
calibration purposes. Cross-talk between adjacent fibres is 
estimated to be less than 5% when using a pure aperture 
extraction (Sanchez 2006). Adjacent fibres in the pseudo- 
slit may come from very different locations on the spatial 
plane (Kelz et al. 2006), minimizing the effect of the cross- 
talk even more (although it does introduce an incoherent 
contamination not important for the present study). 

The V3Q0 grating was used for all the observations, 
covering the wavelength range ~ 3700-7100 A, with a spec- 
tral resolution of FWHM ~ 8 A. Due to the large size of 
NGC 628 (~ 10.5 x 9.5 arcsec) compared to the FOV of the 
instrument a mosaicking scheme was adopted, following the 
experiment by Sanchez et al. (2007b) using the same instru- 
ment. The initial pointing was centered on the centre of the 
galax:y. Consecutive pointings followed an hexagonal pat- 
tern, adjusted to the shape of the PPAK science bundle. 
Each pointing centre is 60 arcsec distant from the previous 
pointing centre. Due to the shape of the PPAK bundle and 
by construction of the mosaic 11 spectra of each pointing. 



28/10/06 
10/12/07 
11/12/07 
12/12/07 
09/08/08 
30/10/08 



1,2,3,9,10,22,23 

1,4,5,6,7,8 

14,15,16,17,18,19,30 

11,12,13,25,29,31,32 

11 

24,26,27,28,29,33,34 



corresponding to one edge of the hexagon, overlap with the 
same number of spectra from the previous pointing. This 
pattern was selected to maximise the covered area, minimise 
large gaps, and allow enough overlap of fibres to calibrate 
exposures taken under different atmospheric conditions. 
For the central pointing we adopted a dithering scheme 
with three positions having offsets (0",0"), (+0. 78", -f 1.68"), 
and (-hO. 78", -1.68") following Sanchez et al. (2007c). This 
dither pattern allows us to cover the gaps between fibres 
and to increase the spatial resolution in this region (where 
more structure is expected) . For each non-dithered pointing 
we obtained 3 exposures of 600s each, and for the central, 
dithered, pointings we obtained 2 exposures of 600s each. 

The observations spanned six nights, distributed over 
three different years. In Fig. 1 the field-of-view of each point- 
ing, labelled with its corresponding identification number, is 
overlaid on top of a _B-band Digital Sky Survey^ Note that a 
substantial fraction of the galaxy is covered by the IFS mo- 
saic. Table 1 gives a log of the observations, including the 
date of observation and the pointings observed each night, 
following the identification numbers shown in Fig. 1. A con- 
centric observing sequence was adopted, starting from the 
inner to the outer regions. The atmospheric conditions var- 
ied between the different observing runs, although in general 
they were clear but non-photometric. The seeing varied be- 
tween ~ 1 and ~ 1.8 arcsec and the median of the seeing 
over all nights was ~ 1.4 arcsec. Different spectrophotomet- 
ric stars were observed during the observing runs, with at 
least two stars observed each night, in order to perform flux 
calibration. In addition to the science pointings, sky expo- 
sures of 300s were taken each night in order to perform a 
proper subtraction of the sky contribution. Additional de- 
tails on the observing strategy can be found in Paper I. 



4 DATA REDUCTION 

Data reduction was performed using R3D (Sanchez 2006), 
in combination with IRAF^ packages and E3D (Sanchez 



^ The Digitized Sky Survey was produced at the Space Telescope 
Science Institute under U.S. Government grant NAG W-2166. 
The images of these surveys are based on photographic data ob- 
tained using the Oschin Schmidt Telescope on Palomar Mountain 
and the UK Schmidt Telescope. The plates were processed into 
the present compressed digital form with the permission of these 
institutions. 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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2004). The reduction consists of the standard steps for fibre- 
based integral-field spectroscopy. A master bias frame was 
created by averaging all the bias frames observed during the 
night and subtracted from the science frames. The different 
exposures taken at the same position on the sky were then 
combined, clipping the cosmic rays, using IRAF routines. 
Then the locations of the spectra on the CCD were deter- 
mined using a continuum illuminated exposure taken before 
the science exposures. Each spectrum was then extracted 
from the science frames. In order to reduce the effects of the 
cross-talk we did not perform a simple aperture extraction, 
which would consist of co-adding the fiux within a certain 
number of pixels of location derived from the continuum il- 
luminated exposure. Instead, we adopted a modified version 
of the Gaussian-suppression technique described in Sanchez 
(2006). 

The new technique assumes a Gaussian profile for the 
projection of each fibre spectrum along the cross-dispersion 
axis. It basically performs a Gaussian fitting to each of the fi- 
bres after subtracting the contribution of the adjacent fibres 
in an iterative process. First, a simple aperture extraction is 
performed with a 5 pixel aperture. This initial guess of the 
flux corresponding to each spectrum is then used to model 
the profiles with a Gaussian function, adopting as centroid 
the location of the peak intensity determined from the con- 
tinuum exposure. The width of the model Gaussian is taken 
from the average width of all the fibre profiles (cr ~ 2 pixels). 
In this first iteration the aperture extracted flux is used as 
the integrated flux of the Gaussian function. This modelled 
profile is then used, for each spectrum, to remove the contri- 
bution of the four adjacent fibre spectra. The resulting clean 
profile is then fitted with a Gaussian function, with the cen- 
troid and width parameters fixed, in order to derive a better 
estimation of the integrated flux. This new flux is used as the 
input for the next iteration of the process. It was found that 
with three iterations the procedure converges, increasing the 
signal-to-noise ratio of the extracted spectra and reducing 
the effect of cross-talk. 

The extracted flux, for each pixel in the dispersion direc- 
tion, is stored in a row-stacked-spectrum flle RSS (Sanchez 
2004). Wavelength calibration was performed using HeHgCd 
lamp exposures obtained before and after each pointing, 
yielding an accuracy of rms ~ 0.3 A. Differences in the rel- 
ative fibre-to-fibre transmission throughput were corrected 
by comparing the wavelength-calibrated RSS science frames 
with the corresponding frames derived from sky exposures 
taken during the twilight. 

Although PPAK is equipped with fibres to sample the 
sky, in most of our pointings the sky fibres are located within 
an area containing significant signal from the galaxy. Accu- 
rate sky subtraction was therefore an important issue in the 
data reduction. For the different nights and pointings we 
adopted different sky subtraction schemes depending on the 
position of the pointing. In some cases, a sufficient number 
of sky-fibres are located in regions free from galaxy emission, 
and it was possible to perform an accurate sky subtraction 
of the individual pointings using these sky spectra. Once a 
certain frame is properly sky-subtracted, the sky spectrum 
of any adjacent frame, observed in the same night and under 
similar weather conditions, can be easily estimated. To do so 
the spectra of the 11 fibres in the sky subtracted frame that 
overlap an adjacent frame (that has not been sky subtracted) 



are combined. This combined spectrum is then subtracted 
from the corresponding combined spectrum of the 11 fibres 
of the adjacent frame. This produces an estimation of the 
sky spectrum in the adjacent frame that is then subtracted 
from all the spectra in that frame. Prior to this subtrac- 
tion it is necessary to visually check that no residual of the 
galaxy is kept in the estimated sky spectrum, which can oc- 
cur if the atmospheric transmission changed substantially 
between the observations of the two frames. 

In other frames the sky observations taken during the 
night could be used to subtract the sky in the science ob- 
servations, especially when the sky and target exposures are 
taken within a few minutes of each other. When the expo- 
sures were more widely separated in time it was necessary to 
combine different sky frames with different weights to derive 
good results. The criterion adopted to decide when a sub- 
traction was good or not was to minimise the residuals in 
the typical emission features of the night spectrum. A thor- 
ough analysis of the sky subtraction and the residuals found 
in the reduction of this galaxy can be found in Paper I. 

After reducing each individual pointing we built a single 
RSS file for the mosaic following an iterative procedure. The 
spectra of each pointing were scaled to those of the previous 
pointing by the average ratio in the 11 overlapping spec- 
tra. Those overlapping spectra were then replaced by the 
average between the previous pointing and the new rescaled 
spectra. The resulting spectra were incorporated into the fi- 
nal RSS file, updating the corresponding position table. By 
adopting this procedure the differences in the spectropho- 
tometric calibration night-to-night and frame-to-frame are 
normalized to that of the first frame used in the process. 
For this reason the mosaic was constructed starting from the 
central frame observed in the night of the 10th of Decem- 
ber 2007, under nearly photometric conditions according to 
the Calar Alto extinction monitor. The final mosaic dataset 
comprises 11094 non-overlapping individual spectra, cover- 
ing a field-of-view of ~6'x7', i.e. ~ 70% of the optical size 
of the galaxy (defined by D25 mag arcsec"^ radius in the 
B-band), and therefore represents the largest spectroscopic 
survey of a single galaxy. 

4.1 Accuracy of the flux calibration 

Although particular care has been taken to achieve the 
best spectrophotometric calibration, there many effects that 
can strongly affect it. Among them, the most obvious are 
the photon-noise from low surface brightness regions of 
the galaxy, the sky-background noise or variations in the 
weather conditions between the time when the spectropho- 
tometric standards and the object were observed. This latter 
effect was reduced by the adopted mosaicking procedure in 
the data reduction, since the photometric calibration was 
renormalized to that of a particularly good night. Less ob- 
vious is the effect of inaccuracies in the sky subtraction. 
However, for low surface brightness regions this is one of the 
most important effects. 

Since it is our goal to provide accurate spectrophoto- 
metric data, we performed a fiux re-calibration on the data. 
To do so we used the flux calibrated broad-band optical im- 
ages from the SINGS legacy survey (Kennicutt et al. 2003). 
In particular, we compared our dataset with the B, V, R and 
Ha images, since they are mostly covered by the wavelength 
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Figure 2. Left: SINGS S-band image of NGC628 used to perform the absolute flux re-calibration. The position of each individual fibre 
of the IFS mosaic is shown overlaid to the real scale. Right: Distribution of the ratio between the flux in the apertures corresponding to 
the PPAK fibres on each of the flux-calibrated broad-band images obtained from the SINGS ancillary data and the corresponding flux 
extracted from the IFU data. Each panel shows the results for the different imaging data: _B-band (top-left panel), V-band (top-right 
panel), -R-band (bottom-left) and Ha (bottom-right panel). The horizontal dashcd-line indicate the one-to-one ratio, which the vertical 
dotted-line indicate the adopted intensity cut for acceptable flux calibration. 



range of our spectra. The photometric calibration of those 
images is claimed to be ~ 5% for the broad-band images 
and ~ 10% for the narrow-band images. They reach a depth 
of ~ 25 mag/arcsec with a signal-to-noise of ~ 10a. There- 
fore, for the structures included in the FOV of our IFS data, 
the photometric errors of the imaging are dominated by the 
accuracy of the calibration, and not by the photon-noise. 

The IFS mosaic was registered to the SINGS optical 
data by matching the coordinates of the galaxy bulge and 
foreground stars in the FOV (see the left-panel of Fig. 2). 
After this process the error in the astrometry of the IFS data 
was estimated to be ~ 0.3" , based on the rms of the differ- 
ences in the centroid of the stars and galaxy bulge. Once 
registered, the mosaic position table was used to extract 
aperture photometry from each broad-band image at the lo- 
cation of each fibre and with an aperture similar to that 
of the fibres. The photometry was transformed to flux (in 
cgs units) by using the counts-to-magnitude prescription in 
the SINGS documentation^, and the zero-points included in 
Fukugita et al. (1995). On the other hand, each spectrum in 
the mosaic was convolved with the corresponding transmis- 
sion curve of the filters indicated before, in order to extract 
a similar flux flbre-to-flbre, based on the IFS data. This pro- 
cedure provides us with 11094 photometric points per band 
to compare between the two datasets. 

The right-panel of Fig. 2 shows the ratio between the 
two sets of photometric points versus the fiux extracted from 
the IFU data, for each of the considered filters. The fig- 
ure shows the typical pattern obtained when comparing the 
fiux ratio between two datasets with different depth, with 
the broad-band images clearly deeper than the spectroscopic 
data, as expected. Down to 0.2 lO^'^'^ erg s~^ cm~^(~5772 
fibres) the ratio is ~ 0.9 for all the fibres, with a standard 



deviation of ~ 0.3 dex. Small differences in the transmission 
curves of the filters used in this calculation and the ones 
used by Fukugita et al. (1995) as well as the astrometric er- 
rors described before introduce uncertainties in the derived 
fiuxes which increase the standard deviation. The uncertain- 
ties from the former source are difficult to estimate. However 
the uncertainties introduced by astrometric errors are esti- 
mated to be at least ~ 10%, by simulating different mosaic 
patterns, moving the location of fibres randomly within 0.3 
arcsec (the uncertainty of our astrometry) and comparing 
the extracted photometry. Based on all these results we es- 
timate our spectrophotometric accuracy to be better than 
~0.2 mag when we apply the re-calibration derived by this 
fiux ratio analysis. The overall re-normalization factor ap- 
plied to the IFS mosaic was 1.15, derived as the average 
of the ratios found on each band. This difference in the 
zero-point of the fiux calibration lies within the expectations 
based on the estimated accuracy of our spectrophotometric 
calibration. 

Fig. 3 shows in the left-panel, the comparison between 
the drift-scan spectrum of the central region of NGG 628 
(blue line), published by Kennicutt et al. (2003)^, and the 
integrated spectrum extracted from the spectrophotomet- 
rically re-calibrated PINGS mosaic (black line), after co- 
adding the spectra within a simulated aperture with the 
same size, location and PA as the SINGS drift-scan. The 
SINGS drift-scan corresponds to a 20" aperture and 70° PA. 
The right-panel of Fig. 3 shows a 100 A width narrow-band 
map of the mosaic of NGC 628 centered at Ha, the red box 
in the centre corresponds to the simulated aperture from 
which the IFS spectrum was extracted. The coordinates, 
size and PA of the aperture were obtained from the header 
of the SINGS data file. Some gaps in between the edges of 
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Figure 3. Left: Comparison between the spectrum extracted for the central region of NGC 628 (shown in the right-panel) using the 
PPAK data (black-line) after co-adding the spectra within a simulated aperture with the same size, location and PA as the SINGS 
legacy survey drift-scan (blue-line, in the online version of the paper). The bottom-panel shows the ratio between the two spectra. Right: 
Narrow-band map of the PINGS mosaic of NGC 628 centered at Ho, the red box in the centre corresponds to the simulated aperture 
from which the IFS spectrum was extracted. Despite the different techniques used to derive both spectra, there is a clear agreement 
between them in most of the wavelength range. 



the pointings are seen in this map, they are due to the re- 
centering of the individual pointings after comparing with 
the broad-band images as discussed in above. As expected 
from the spectrophotometric re-calibration of the mosaic, 
there is a very good agreement between both datasets; both 
the general shape of the spectra and the strength of the 
spectral features match well. The left-panel of Fig. 3 shows 
in the lower part the relative difference between the two 
spectra, which is consistent with a null difference (in aver- 
age) and within the error of the absolute flux calibration 
for most of the spectral range. As simulations with different 
position, apertures and mosaic versions showed, the small 
deviations in the continuum are due to the presence of a 
foreground star within the FOV of the simulated slit. The 
small disagreement at wavelengths shorter than 4000 A is 
expected due to the degradation of the CCD sensitivity in 
the blue, as explained in detailed in Paper I, where further 
comparisons with previous spectrophotometric data can be 
found. A sligth missmatch of the wavelength resolution at 
the edges of the spectra is also noticed, being in the range 
of the expectations for such comparisons. 

Fig. 4 shows the comparison between two reconstructed 
V-band images of NGC 628. The image on the left was cre- 
ated after interpolating the aperture photometry extraction 
of the SINGS broad-band image. The image on the right 
was derived from the flux re-calibrated IFS dataset, once 
each spectrum corresponding to a particular fibre was con- 
volved with the filter response, as described before. In or- 
der to create a regularly gridded image, the data were in- 
terpolated using E3D (Sanchez 2004), adopting a natural- 
neighbor non-linear interpolation scheme, and a final pixel 
scale of 1" /pixel. The areas not fulfilling our criteria for 
accurate spectrophotometry (Fv <0.210~'^^ erg s"'^ cm^^) 
were masked. 



5 ANALYSIS AND RESULTS 

In order to extract physical properties of the galaxy from the 
dataset it is necessary to perform different kinds of analy- 
ses. In particular, for each spectrum we need to identify the 
emission lines of ionized gas, and decouple this emission from 
the underlying stellar population. Particular care has to be 
taken in this decoupling technique, since some of the emis- 
sion lines (e.g. H/3) may be strongly affected by underlying 
absorption features. 

The decoupling of the stellar population from the emis- 
sion lines was performed adopting a scheme summarized 
here: (i) A set of detected emission lines was identified in the 
integrated spectra of the stronger H II regions in the outer 
part of the galaxy, (ii) For each spectrum in the dataset, the 
underlying stellar population was fitted by a linear combina- 
tion of a grid of single-stellar population (SSP) templates, 
after correcting for the appropriate systemic velocity and 
velocity dispersion (including the instrumental dispersion 
which dominates the total observed dispersion), and tak- 
ing into account the effects of dust attenuation. A spectral 
region of 30 A width around each detected emission line was 
masked prior to the linear fitting, including also the regions 
around the sky-lines (Sanchez et al. 2007a). (iii) Once we 
derived a first approximation of the spectrum of the under- 
lying stellar population, this was subtracted from the origi- 
nal spectrum to obtain a pure emission-line spectrum, (iv) 
To derive the intensity of each detected emission line, each 
of these clean spectra in the dataset was fitted to a single 
Gaussian function per emission line plus a low order poly- 
nomial function. Instead of fitting emission lines over the 
entire wavelength range simultaneously, for each spectrum 
we extracted shorter wavelength ranges that sampled one or 
a few of the analysed emission lines, in order to characterise 



8 S. F. Sanchez et al. 



-100 



-200 



SINGS 



200 



100 - 



b 

CD 




h 



-100 



-200 



Aa (arcsec) 



PINGS 



CO 200 



100 - 




Aa (arcsec) 



Figure 4. Reconstructed y-band images of NGC 628. Left: image created after interpolating the aperture photometry extraction of the 
SINGS broad-band image. Right: interpolated image derived after multiplying each spectrum of the IFS mosaic with the filter response 
curve. Both images were created with a regular grid of 1" /pixel. The areas not satisfying our criteria of accurate spectrophotometry were 
masked. The perpendicular lines are centered at the mosaic's reference point. Offsets are in arcsec. 



the residual continuum with a simple polynomial function, 
and to simplify the fitting procedure. When more than one 
emission line was fitted simultaneously, their systemic veloc- 
ities and FWHMs were forced to be equal (since the FWHM 
is dominated by the instrument resolution), in order to de- 
crease the number of free parameters and increase the accu- 
racy of the deblending process (when required) . (v) Finally, 
for each spectrum in the dataset a pure gas-emission spec- 
trum was created, based on the results of the last fitting pro- 
cedure, using only the combination of Gaussian functions. 
This model was then subtracted from the original data, spec- 
trum by spectrum, to produce a dataset of gas-free spec- 
tra. These spectra are then fitted again by a combination 
of SSP, as described before (but without masking the spec- 
tral range around the emission lines, in this case), deriving 
the lummosity-weighted age, metallicity and dust content of 
the composite stellar population. The Appendix A gives a 
detailed description of the fitting procedure, indicating the 
basic algorithms adopted, and including estimates of the ac- 
curacy of the SSP-based modeling and the derived parame- 
ters based on simulations. The 2D maps of the emission line 
intensities and physical properties shown in the following 
sections were constructed based on the pure gas-emission 
mosaics described above. 



5.1 Integrated spectrum 

A particularly interesting use of IFS datasets is the com- 
bination of the observed spectra to produce an integrated 
spectrum of the object, using the IFU as large aperture spec- 
trograph. This high signal-to-noise integrated spectrum can 
be used to derive, for the first time, the real average spectro- 
scopic properties of the galaxy, as opposed to previous stud- 



ies that attempted to describe its average properties by the 
analysis of individual spectra taken at different regions. The 
most similar approach would be the spectrum derived by us- 
ing a drift-scanning technique (e.g. Moustakas & Kennicutt 
(2006b) and part of the ancillary data of the SINGS sur- 
vey), although in those studies (specially for the latter) the 
fraction of galaxy covered by the spectra was much less than 
the for the spectrum presented here. Another advantage of 
the use of an IFU with respect to the drift-scan technique is 
that the former allows a comparison between the integrated 
and the spatially resolved properties of the galaxy. 

A good amount of the fibres in the IFS mosaic of 
NGC 628 do not contain enough signal-to-noise or do not 
contain signal at all (i.e. spectra with a flat continuum con- 
sistent with zero-fiux), as the fibres were sampling regions 
where the intrinsic fiux from the galaxy is low or null (e.g. 
borders of the mosaic, intra-arms regions, etc). In order to 
get rid of spectra where no information could be derived, we 
obtained a clean version of the IFS mosaic of NGC 628 by 
applying a fiux threshold cut choosing only those fibres with 
an average fiux along the whole spectral range greater than 
10~^® ergs~^ cm~^ A~^. Furthermore, bad fibres (due to 
cosmic rays and CCD cosmetic defects) and foreground stars 
(10 within the observed field-of-view of NGC 628) where re- 
moved from the mosaic. This procedure resulted in a refined 
mosaic version which includes only those regions with high- 
quality spectrophotometric calibration. The total number of 
spectra in the clean IFS mosaic version of NGC 628 is 6949. 
The top-panel of Fig. 5 shows the integrated spectrum of 
NGC 628 derived by co-adding the spectra within the clean 
IFS mosaic version of the galaxy. 

The integrated spectrum of NGC 628 shows a character- 
istic stellar continuum with absorption features and emission 
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Figure 5. Integrated spectrum of NGC 628. The red solid 
squares in the top-panel indicate the integrated flux derived from 
the B, V, R and Ho images. The middle-panel shows the SSP 
model fitting (red-line, in the online version) to the spectrum, 
the light-blue bands correspond to the spectral regions masked 
during the fitting. The bottom panel shows the residual after 
subtracting the model from the original spectrum, the detected 
emission lines have been labeled. 



lines superimposed. Ha, H/3, [O II] A3727, [O III] A5007, the 
[Nil] AA6548,84 and AA6717,31 doublets are clearly identi- 
fied. Less obvious are the H7 and [O III] A4959 lines. Sky 
residuals are also present in the spectrum, especially the 
[O I] A5577 and [Na I] A5893 lines. The red-solid squares (in 
the online version) correspond to the integrated flux derived 
from the B, V, R and Ha images of the SINGS ancillary 
data, obtained during the spectrophotometric re-calibration, 
using the same apertures of the IFU data and co-adding 
them in a similar way as the integrated spectrum. The posi- 
tion of these data-points with respect to the continuum (in 
the case of the broad-band images) and the peak of the Ha 
line in the spectrum, corroborates the accuracy of the abso- 
lute flux calibration. The co-added region comprises ~90% 
of the total flux of the galaxy in 1/-band, as estimated from 
corresponding SINGS image. An analysis similar to that de- 



scribed above was performed on this integrated spectrum to 
derive the main physical properties of both the ionized gas 
and the stellar population. 



5.1.1 Integrated stellar populations 

The middle-panel of Fig. 5 shows the best SSP model fit 
to the integrated spectrum superimposed in red colour. The 
light-blue bands correspond to the spectral regions masked 
during the fitting as explained above. They coincide with 
the position of the strongest redshifted emission lines and 
regions of bright sky residuals. Note the strength of the un- 
derlying stellar absorptions in the Balmer lines. The multi- 
component SSP model matches accurately the continuum of 
the integrated spectrum, within an error of the ~3%. Most 
of the discrepancies are in regions clearly dominated by im- 
perfections in the sky-subtraction, due to the strength of the 
night sky lines (e.g. at ~ A5577, see Sanchez et al. (2007a) 
and Paper I). In order to evaluate the possible effects of 
these residuals in the derivation of the main properties of the 
stellar population the analysis was repeated for the wave- 
length range between A4100 and A5400 A, where there are 
no such strong atmospheric features. For this wavelength 
range, the model matches the continuum within a ~2%. In 
the UV-regime the errors are sometimes larger than this 
value. 

Table 2 lists the luminosity-weighted age, metallicity 
and dust attenuation of the best fitting model for the two 
cases considered here: one using the entire spectral range in 
the fitting process and the other using the reduced range de- 
scribed before. There is a very good agreement between the 
results derived in both cases. Based on the results from the 
simulations (see Appendix), we expect an error of 10-20% 
in the derivation of the Age, and of a ~5% in the derivation 
of the metallicity. These errors do not take into account the 
systematics, due to the applied algorithm, the current code 
and the templates adopted, or the known degenerancies in 
the derivation of the age and the metallicity, and therefore 
the expected discrepancies with previous published results 
are much larger than the ones estimated on the basis of 
the simulations. From these fits we may conclude that the 
average stellar population of NGC 628 is dominated by an 
old component of ~ 9 Gyrs, but with a subsolar metallicity 
{[Z/H] ~ —0.45). Although synthesis modeling is nowadays 
widely used, we have to take into account the age- metallicity 
degeneracy problem that plagues most spectral fitting tech- 
niques. The most frequently used technique to derive the 
age and metallicity of the stellar population in galaxies is to 
measure certain line strength indices, such as the Lick/IDS 
index system (e.g. Burstein et al. 1984; Faber et al. 1985; 
Burstein et al. 1986; Gorgas et al. 1993; Worthey 1994). For 
this purpose, one generally tries to use a combination of in- 
dices that are most orthogonal in the parameter space (i.e. 
age and metallicity). In order to cross-check the results based 
on the fitting procedure, we measured the equivalent width 
of H/3 (primarily sensitive to the age) and Mgfe (primarily 
sensitive to the metallicity), of the integrated spectrum. 

The equivalent widths of the absorption lines were de- 
rived using the bandpass definitions from the Lick index sys- 
tem revised by Trager et al. (1998), shifted to the redshift 
of the object, as described in Sanchez et al. (2007c). To de- 
rive the age and metallicity from the measured indices, we 
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Table 2. Properties of the average stellar population. 



Analysis method 


Age 

(Gyr) 


Metallicity 
Z [Z/H] 


(mag) 


SSP fit, 3700-6800 A 
SSP fit, 4100-5400 A 
H/3 vs. Mgb indices 


8.95 
8.40 
9.78 


0.007 
0.007 
0.008 


-0.44 
-0.44 
-0.42 


0.4 
0.8 



* Dust continuum attenuation. 

have adopted the model grid from Thomas et al. (2003), im- 
plemented in the rmodel code^. The resulting estimates of 
the age and metalhcity, based on the absorption line index 
analysis, are listed in Table 2. Despite the strong conceptual 
differences between these methods and the fitting technique 
described before, the results are very consistent with the val- 
ues obtained in both the full and reduced wavelength range 
SSP fitting. 

5.1.2 Integrated properties of the ionized gas 

The bottom-panel of Fig. 5 shows the pure emission line 
spectrum of NGC 628 obtained after subtracting the under- 
lying stellar population from the integrated spectrum. As 
expected, the spectrum is dominated by a set of emission 
lines, plus a residual continuum consistent with a zero av- 
erage intensity. This procedure reveals additional emission 
lines, like H(5 and He. All the detected lines have been la- 
beled with their standard notation. Each of these emission 
lines was fitted with a single Gaussian function, as described 
before, in order to derive their strengths. Fig. 6 shows an ex- 
ample of the fits in the wavelength range between 6500 A and 
6800 A. The black solid-line corresponds to the pure emis- 
sion line spectrum of NGC 628, the red dashed-line shows 
the best fitted model, consisting of six Gaussian functions, 
assuming a single Gaussian fit to each of the emission lines 
detected in this wavelength range: the [Nil] doublet. Ha, 
He I A6678 (very faint) and the [S II] doublet. In this partic- 
ular case both the systemic velocity and velocity dispersion 
were forced to be the same for all the emission lines, in order 
to increase the accuracy of the derived parameters. The ratio 
between the two [Nil] lines included in the spectral range 
were fixed to the theoretical value (Osterbrock & Ferland 
2006). A similar procedure was applied to the rest of the 
emission lines. By adopting this procedure, it is possible to 
accurately deblend the different emission lines. 

Table 3 lists the emission line intensities derived us- 
ing this procedure, including for each of the detected emis- 
sion lines its standard identification, the laboratory rest- 
frame wavelength, and the estimated intensity normalised 
to the observed flux of H/? in units of 10~^^ ergs~'^cm~^. 
They are shown in the columns labeled as F(A)/H/3. 
The observed [O III] ratio shown in Table 3 is consistent 
with the well-known theoretical value between these lines 
(Storey & Zeippen 2000). The associated la errors are solely 
due to the statistical uncertainty (Tstat in the measurement 
of the fiux intensity. Thereafter, the observed line intensi- 
ties were corrected for reddening using the Balmer decre- 
ment according to the reddening function of Cardelli et al. 




6600 



6700 



Wavelength (A) 

Figure 6. Example of the fitting procedure applied to derive 
the intensity flux of the detected emission lines. The black solid- 
line corresponds to the pure emission line spectrum of NGC 628, 
covering the wavelength range between 6500 and 6800 A. The 
red dashed-line shows the best fitting model describing the emis- 
sion lines, which comprises a single Gaussian function for ea.ch of 
them. The green dotted-line shows the residual spectrum, once 
the previous model has been subtracted. 



Table 3. Integrated line intensities for NGC 628. The first col- 
umn corresponds to the emission line identification, with the rest- 
frame wavelength, the second one corresponds to the adopted 
reddening curve normalized to H/3. The F(A)/H/3 column corre- 
sponds to the observed flux, while the I(A)/H/9 to the reddening 
corrected values; normalised to H/3. The values in parenthesis cor- 
respond to the la errors calculated as explained in the text. The 
observed flux in H/3 is in units of 10~^^ ergs~^cm~^. The last 
row shows the number of fibres from which the integrated spectra 
was extracted. 



Line 


/(A) 


F(A)/H/3 


I(A)/H/3 


[O II] A3727 


0.32 


1.494 


0.046) 


2.138 


0.320) 


H8 + He I A3889 


0.29 


0.158 


0.032) 


0.217 


0.054) 


H(5 A4101 


0.23 


0.165 


0.032) 


0.214 


0.051) 


H7 A4340 


0.16 


0.511 


0.023) 


0.609 


0.083) 


H/3 A4861 


0.00 


1.000 


0.014) 


1.000 


0.052) 


[O HI] A4959 


-0.03 


0.111 


0.010) 


0.108 


0.017) 


[O HI] A5007 


-0.04 


0.329 


0.010) 


0.315 


0.040) 


He I A5876 


-0.20 


0.299 


0.012) 


0.239 


0.031) 


[N II] A6548 


-0.30 


0.383 


0.013) 


0.276 


0.035) 


Ha A6563 


-0.30 


3.998 


0.041) 


2.870 


0.352) 


[N II] A6584 


-0.30 


1.111 


0.016) 


0.795 


0.098) 


[SII] A6717 


-0.32 


0.619 


0.011) 


0.435 


0.054) 


[S II] A6731 


-0.32 


0.393 


0.010) 


0.276 


0.035) 


[O HI] A5007/A4959 




2.96 


(0.28) 


2.92 


(0.58) 


[S II] A6717/A6731 




1.57 


(0.14) 


1.58 


(0.28) 


F(H/3) A4861 






1549.4 






c(H/3) 




0.48 


(0.05) 






Av 






1.04 






Extraction fibres 






6949 







N. Cardiel, http://tinyurl.com/rmodel 



Wide-field IFS of NGC 



11 



(1989), assuming R = Av/E{B - V) = 3.1. Theoretical 
values for the intrinsic Balmer line ratios were taken from 
Osterbrock & Ferland (2006), assuming case B recombina- 
tion (optically thick in all the Lyman lines), an electron den- 
sity of Ue = 100 cm~^ and an electron temperature Te = lO** 
K. The last column of Table 3 shows the reddening-corrected 
emission line fluxes for the integrated spectrum, designated 
by I(A)/H/3. These flux ratios can be used to derive the aver- 
age properties of the ionized gas in the galaxy. The adopted 
reddening curve normalized to H/3, /(A), is shown in the 
second column of the table. Although higher-order Balmer 
lines were detected in the integrated spectrum, no Balmer 
lines beyond H7 were used for the determination of c(H/3), 
as the associated error of the measurement of these lines in 
the residual spectrum yielded high uncertainties in the com- 
puted c(H/3), due to their low signal-to-noise; therefore the 
c(H/3) value was derived using the Hq/H/? ratio solely. No 
auroral lines were detected either in the integrated or resid- 
ual spectrum. The low strength of the [O III] lines suggests 
a relatively high metallicity for the integrated abundance of 
NGC 628, as it will be shown hereafter. 

Although particular care has been taken in the flux cal- 
ibration of the spectra within the mosaic, the absolute flux 
intensity listed in Table 3 has to be treated with caution. 
On one hand, not all of the galaxy surface has been cov- 
ered by our IFS observations, as can be seen in Fig. 1. In 
particular some of the brighter H II regions, located to the 
east of NGC 628 are not included in the FOV of our mo- 
saic. This causes to underestimate the intensity of all the 
lines. By a rough estimation, based on the D25 optical radius 
(_B-band), we consider that our IFS Mosaic covers ~ 70% 
of the galaxy size. On the other hand, the central fibre- 
bundle of PPAK has a filling factor of ~ 65%, as mentioned 
in Section 3, which leads to a corresponding underestimation 
of the integrated flux. In the particular case of this mosaic, 
a dithering scheme necessary to compensate for this incom- 
plete sampling was adopted only for the central pointing, 
where the emission lines are in general weak. All together, 
we consider that it is necessary to apply a correction of a 
factor ~ 2.2 to take into account the aperture and sampling 
effects described before. As mentioned before, NGC 628 has 
been extensively studied before in several publications. In 
particular, different authors have reported on the Ha in- 
tensity flux, using different procedures, from photoelectric 
photometers to narrow-band imaging. Table 4 lists a sum- 
mary of these published values, together with the value de- 
rived from our integrated spectrum and considering the flux 
correction factor described above. Despite the different bi- 
ases introduced by the different methods, there is substantial 
agreement between the previously published results and our 
reported value for the integrated Ha flux of the galaxy. 

The derived dust extinction, Av = 1.04, larger than 
the one derived from the analysis of the stellar populations 
(Av ~ 0.4 mag). This result is not surprising, since both 
methods sample different regions of the galajcy. While the 
underlying continuum is dominated by the stellar compo- 
nents of the central regions, clearly brighter, the ionized gas 
spectrum is dominated by the star-forming regions in the 
spiral arms. These latter regions are known to be more at- 
tenuated by dust, due to star forming process (e.g. Calzetti 
2001). Indeed, the extinction law derived by Calzetti (1997) 
for star-forming galaxies shows that the typical extinction in 



Table 4. Comparison between different H« fluxes reported for 
NGC 628 in tlie literature. Fluxes in units of 10~^^ era; s~^ crn"^. 



Flux 



Reference 



1.07 Kennicutt 1983 

0.87 Young et al. 1996 

1.51 Hoopes et al. 2001 

1.05 Marcum et al. 2001 

1.02 Kennicutt et al. 2008 

1.14 Current study 



the emission lines of these objects is approximately double 
that in their stellar continuum. 

The integrated flux of Ha and [O II] A3727 can be used 
to determine a rough value of the global star formation rate 
(SFR) in this galaxy. The intensities of both lines were cor- 
rected by dust extinction, adopting the Av value and the 
aperture correction mentioned before. Absolute luminosi- 
ties were derived by assuming a standard ACDM cosmology 
with Ho=70.4, Q,rn=0-268 and S1a=0.73, and a luminosity 
distance of 9.3 Mpc (Hendry et al. 2005). The derived lumi- 
nosities are Lhq, ~3.08 and L[oii] '^ 2.30, in units of lO''^ 
erg s~^. The values of the SFR were derived adopting the 
classical relations by Kennicutt (1998), obtaining SFR ~2.4 
and 3.2 Mq yr~^, based on the Ha and [Oil] luminosities 
respectively. 

Different possible mechanisms can be responsible for the 
ionization in emission line galaxies. The nature of the ioniza- 
tion can be derived from ratios of the usual diagnostic lines 
(BPT, Baldwin et al. 1981; Veilleux & Osterbrock 1987). 
Based on the values listed in Table 3, it was found that 

logio([NlI] A6584/Ha) 0.56 and logio([0lll] A5007/H/3) 

~ —0.5. These line ratios correspond to the expected 
values for star forming galaxies and/or HII regions (e.g. 
Sanchez et al. 2005, 2007c), far away from any of the bound- 
ary regions in the [O III] /H/3 vs. [N II] /Ha diagnostic di- 
agram. Therefore, it is clear that the dominant ionization 
mechanism in the integrated spectrum of NGC 628 is due 
to hot OB stars (HII regions), as expected, since no previ- 
ous study has reported any kind of nuclear activity in this 
galaxy. 

The (volume-averaged) ionization parameter, defined as 
the ratio of the density of ionizing photons to the parti- 
cle density; u = S^' , where Qgo is the flux of ionizing 
photons produced by the exciting stars above the Lyman 
limit, Rs is the radius of the Strmgren sphere, n the number 
density of hydrogen atoms, and c the speed of light. This 
parameter determines the degree of ionization at any par- 
ticular location within the nebula. As discussed in the next 
section, many of the empirical methods commonly used to 
derive the chemical abundance of a star-forming region are 
sensitive to this parameter, and for some ranges of metallic- 
ity, they are not useful unless the ionization parameter can 
be constrained within a small range of possible values. 

The ionization parameter is best determined using the 
ratios of emission lines of different ionization stages of the 
same element. In general, the larger the difference in ioniza- 
tion potentials of the two stages, the better the constrain. A 
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commonly used ionization parameter diagnostic is based on 
the ratio [O II] /[O III] = A3727 / (A4959 + A5007). However, 
this ratio is not only sensitive to the ionization parameter, 
but is also strongly dependent on metallicity. Another possi- 
bility is to use the [S III] A9069 and/or [S III] A9532 together 
with the [SII] A6717, A6731 emission lines. The [S II]/[S III] 
ratio provides a more reliable useful ionization parameter di- 
agnostic (Diaz et al. 1991). Considering the available lines 
within the spectral range of the PINGS observations (Diaz 
1994), a good approximation to the ionization parameter 
can be determined from the [O Il]/H/3 or [S II] /Ha ratio , 



NGC 628 integrated spectrum 



log u = -0.80 log ([O II] / [O III]) - 3.02, 



(1) 



after Diaz et al. (2000). Its uncertainty is estimated to be 
±0.2 dex. The derived logu for integrated spectrum of 
NGC 628 is -3.58, which corresponds to a low value of the 
ionization parameter (--^ 10~*). 

5.1.3 Integrated oxygen abundance 

Galaxies in the local universe have been used as an anchor 
point to determine the evolution of the metallicity (based 
on the gas phase oxygen abundance) , along different cosmo- 
logical periods. The global metallicity of a given galaxy is 
represented by its oxygen abundance. The rest of the ele- 
ments vary within the same fractions as those found in the 
Sun. However, while the determination of this observable 
at high redshift normally describes the average value of the 
galaxies as a whole (due to aperture effects), in the local uni- 
verse most of these determinations are based on studies of a 
number of discrete HII regions. The integrated spectrum of 
NGC 628 can be used to perform an integrated abundance 
analysis in a consistent way as the studies performed over 
high redshift galaxies, with the advantage than the results 
of the integrated study can be compared with the abun- 
dances of resolved regions within the galaxy. In this section, 
we present a chemical abundance analysis of the integrated 
spectrum of NGC 628 based on a suite of different empirical 
abundance diagnostic methods. A more detailed comparison 
between the integrated and spatially resolved abundances 
will be presented in a subsequent paper (Rosales-Ortega et 
al., in preparation). 

Given the limitations imposed by the non-detection 
of auroral temperature sensitive emission lines (such as 
[O III] A4363 or [N II] A5755) in HII regions of low-excitation 
and/or low surface brightness, empirical methods based on 
the use of strong, easily observable optical lines have been 
developed throughout the years. Although abundances de- 
rived in this way are recognised to suffer considerable un- 
certainties, they are believed to be able to trace large-scale 
trends in galaxies. Following these ideas, several abundance 
calibrators have been proposed involving different emission- 
line ratios and have been applied to determine oxygen abun- 
dances in objects as different as individual HII regions 
in spiral galaxies, dwarf irregular galaxies, nuclear star- 
burst, and the integrated abundances of spiral galaxies (e.g. 
Kobulnicky et al. 1999; Pilyugin et al. 2004, and references 
therein) . 

In order to derive the characteristic chemical abun- 
dance of NGC 628, a set of empirical calibrators were ap- 
plied to the integrated spectrum. Different abundance es- 
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Figure 7. Comparison of the integrated oxygen abundance of 
NGC 628 derived for different estimators. Red-solid points corre- 
spond to iJ23 calibrators, green-points to index-empirical meth- 
ods, and blue-points to the different methods proposed by Pilyu- 
gin and collaborators. The open circles correspond to an arbitrary 
-0.5 dex offset of the R23 based methods. The horizontal dashed- 
line correspond to the oxygen solar value. 



timators were chosen in order explore the effect of a par- 
ticular calibration depending on the physical properties of 
the galaxy. By far, the most commonly used such ratio is 
([Oil] A3727 + [QUI] A4959, A5007)/H/?, known as the i?23 
method (Pagel et al. 1979). The logic for the use of this ra- 
tio is that it is not affected by differences in relative ele- 
mental abundances, and remains essentially constant within 
a given giant HII region despite variations in excitation 
(Diaz et al. 1987). However, there is an ambiguity inherent 
to this method since there are two values of abundance cor- 
responding to a given value of R23, i.e. the lower branch in- 
creases with increasing abundance, while the upper branch 
shows an opposite behavior, i.e. 7^23 decreases with increas- 
ing abundance. Another drawback of using i?23 (and many 
of the other emission-line abundance diagnostics) is that it 
depends also on the ionization parameter. From this cate- 
gory, two methods were applied: the McGaugh 1991 (here- 
after M91) and the Kobulnicky & Kewley 2004 (hereafter 
KK04) calibrations, which are theoretical methods based on 
photoionization models. Both take into account the ioniza- 
tion parameter to produce an estimate of the metallicity. 
The estimated accuracy of these methods is ~ ±0.15 dex. 

In order to discriminate between the different _R23 
branches, we used the [Nil] A6584/[0II] A3727 ratio fol- 
lowing the prescriptions by Kewley & Dopita (2002). The 
[Nil] /[On] ratio is not sensitive to the ionization parameter 
to within ±0.05 dex, and it is a strong function of metallicity 
above log([Nll]/[0 II]) > -1.2, where the division between 
the upper and lower branches occurs. The log([Nll]/[0 II]) 
value for the integrated spectrum of NGC 628 is —0.43, i.e. 
indicating that the R2Z value of the integrated spectrum for 
NGC 628 (log i?23 = 0.41) corresponds to the upper branch 
of the i?23 relation. 

Another subset of estimators was chosen from 
the category of empirical strong-line methods, they 
correspond to the N2 calibration (first proposed by 
Storchi-Bergmann et al. 1994, but using the definition af- 
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Table 5. Integrated oxygen abundances for NGC 628 in units of 12+log(0/H). The columns designations correspond to the 
following abundance calibrators: M91, McGaugh (1991); KK04, Kobulnicky & Kewley (2004); N2, Denicolo et al. (2002); 
03N2, Pettini & Pagel (2004); PT05, Pilyugin & Thuan (2005); (0/H)fi, ff-T^ method (as explained in the text). 



log R23 


M91 


KK04 


N2 


03N2 


PT05 


(0/H)ff 


0.41 


8.90 ± 0.14 


9.06 ± 0.11 


8.62 ± 0.16 


8.71 ± 0.11 


8.32 ± 0.15 


8.55 ± 0.09 



ter Denicolo et al. 2002), and the 03N2 calibration (first 
proposed by AUoin et al. (1979), but using the definition by 
Pettini & Pagel 2004). These two indices have the virtue 
of being single- valued, however, they are affected by the 
low-excitation line [Nil] A6584, which may arise not only 
in bona-fide HII regions, but also in the diffuse ionized 
medium, which is an issue for spectra integrated within ex- 
tended regions, such as the integrated spectrum of NGC 628. 
The estimated uncertainty of the derived metallicities is 
~ ± 0.2 dex. 

The Pilyugin & Thuan 2005 (hereafter PT05) calibra- 
tion is based on an updated version of the Pilyugin (2001) 
estimator, obtained by empirical fits to the relationship be- 
tween R23 and Te metallicities for a sample of HII regions. 
This estimator was also considered to determine the in- 
tegrated abundance of NGC 628. This calibration includes 
and excitation parameter P that takes into account the ef- 
fect of the ionization parameter. The PT05 calibration has 
two parametrizations corresponding to the lower and up- 
per branches of 7?23. As in the case of the M91 and KK04 
calibrators, the [Nil] /[Oil] ratio was used to discriminate 
between the two branches of the 7?23 relation. 

The last strong-line empirical method considered is 
a combination of the fiux-flux (or ff-relation) found by 
Pilyugin (2005) and parametrised by Pilyugin et al. (2006), 
the i2 — is relation between the O"*" and O^"'" zones elec- 
tron temperatures for high-metallicity regions proposed by 
Pilyugin (2007), and an updated version of the Te-based 
method for metallicity determination (Izotov et al. 2006). 
According to these authors, the combination of these meth- 
ods solves the problem of the determination of the electron 
temperatures in high-metallicity HII regions, where faint au- 
roral lines are not detected. However, the abundances de- 
termined through this method rely on the validity of the 
classic Te method, which has been questioned for the high- 
metallicity regime in a number of studies by comparisons 
with HII region photoionization models (Stasiriska 2005). 
The abundances derived through this method will be re- 
ferred as (0/H)fi or ff-T^ abundances. 

The advantages and drawbacks of all these differ- 
ent calibrators have been thoroughly discussed in the 
literature (e.g. Pagel et al. 1980, Kennicutt & Garnett 
1996, Kewley & Dopita 2002, Perez-Montero & Diaz 2005, 
Kewley & Ellison 2008). Regrettably, comparisons among 
the metallicities estimated using these methods reveal large 
discrepancies. They are usually manifested as systematic off- 
sets in metallicity estimates, with high values corresponding 
to theoretical calibrations and lower metallicities estimated 
by electron temperature metallicities, with offsets as large 
as 0.5 dex in log(0/H) units (Liang et al. 2006). 

The calculated oxygen abundances (in units of 
12 -f log(0/H)) derived for each of the calibrators listed 
above are shown in Section 5.1.3. The uncertainties cor- 



respond to the Icr error found by propagating the errors 
through a Monte Carlo simulation by using Gaussian dis- 
tributions with a width equal to the errors of the emission 
line intensities, modulated by recomputing the distribution 
500 times. These abundances are plotted in Fig. 7. The red- 
solid points (in the online version) correspond to the _R23 
methods, the green points to the index-empirical methods, 
and the blue points to the different methods proposed by 
Pilyugin and collaborators. The red-open circles correspond 
to an arbitrary -0.5 dex offset of the _R23 methods (these 
offsets are included given the well-known systematic off- 
set of the i?23 theoretical-based calibrations). The horizon- 
tal dotted-line corresponds to the oxygen solar abundance 
12 + log(0/H)=8.7 (Scott et al. 2009). 

The two -R23 methods derive super-solar metallicity val- 
ues, while the PT05 oxygen value correspond to the lowest 
metallicity obtained. The index-empirical methods, N2 and 
03N2, stand in between the metallicities derived through 
the i?23 and the PT05 -I- ff-T^ methods. Interestingly, the 
_R23 calibrations metallicities which were shifted to a lower 
value (red-open circles), are in close agreement (within the 
errors) to the pure-empirical and the ff-Te based metallici- 
ties. The mean oxygen abundance derived from all the meth- 
ods is 12 + log(0/H)=8.69 ± 0.31. 

Previous studies have determined the oxygen abun- 
dance in this galaxy based on spectroscopic imaging spec- 
trophotometry observations of different H II regions. In par- 
ticular, McCall et al. (1985) analysed a sample of HII re- 
gions within a radius of ~ 200 arcsec, basically coincident 
with the area sampled by our IFS survey. They reported 
a range of oxygen abundances between 12 -I- log(0/H) 
~ 8.7 — 9.3 (by employing a method based on the -R23 index), 
with a considerable decline from the inner to the outer parts. 
Belley & Roy (1992) obtained reddenings, H/3 equivalent 
widths, diagnostic line ratios and metallicities for 130 H II re- 
gions by the implementation of imaging spectrophotome- 
try. They derived an abundance gradient of NGC 628 based 
on the [Qui] /H/3 empirical calibrator (Edmunds & Pagel 
1984). Their values range between 12 + log(0/H) ~ 8.4-9.2, 
covering a large baseline in galactocentric distances (up to 
p ~ 2peff). Although this was not a strict spectroscopic 
study, given the number of H II regions analysed, the work of 
Belley & Roy 1992 stood up to now as the most complete 2D 
description of the emission line chemistry of NGC 628. Sub- 
sequently, van Zee et al. (1998) reported the oxygen abun- 
dances of the HII regions in an outer ring between ~ 150 
arcsec and ~ 300 arcsec radius. They found that the de- 
cline in the abundance continues, and the oxygen abun- 
dance ranges between 12 ~f log(0/H) ~8.10 and 8.95 (using 
a modified version of the M91 calibrator). More recently, 
Castellanos et al. (2002b) observed a reduced set of HII re- 
gions in the optical and near-infrared where they were able 
to measure temperature-sensitive emission lines. They re- 



14 S. F. Sanchez et al. 




O 






o 



CD 
O 




50 



Radius (arcsec) 



100 



Figure 8. Left panel: radial distribution of the luininosity-wcight age (in Gyr), of the stellar population based on the fitting procedure 
using the combination of multiple SSPs. Right panel: similar distribution for the luminosity-weight metallicity, relative to the solar one 
(Z/H). 



ported an average oxygen abundance of 12 + log(0/H) 
~ 8.23. However these HIT regions are located beyond 
the FOV of our IFS mosaic, at galactocentric distances 
where, given the well-known radial metallicity gradient of 
this galaxy (Ferguson et al. 1998), the oxygen abundance 
is expected to be lower than the integrated metallicity de- 
rived at inner radii. Therefore we cannot compare the inte- 
grated 0/H abundance derived in this work with the values 
of Castellanos et al. (2002b) due to the non-coincident ge- 
ometry. On the other hand, the range of oxygen abundances 
reported by the previous spectroscopic and imaging spec- 
trophotometry studies agree perfectly with the integrated 
abundances derived through the R23, methods and with the 
mean integrated abundance of NGC 628. 

The integrated properties of the ionized gas derived in 
this chapter need to be compared with the resolved proper- 
ties in order to analyse the validity of the results obtained 
from the integrated analysis, taking into account different 
effects, such as the extinction or the contribution of the dif- 
fuse interstellar emission. These points will be addressed in 
the following sections. 



5.2 Spatially resolved properties of the galaxy 

In this section we present the results obtained by apply- 
ing the spectra fitting technique described in Section 5 to 
each individual spectrum of the dataset. As indicated in the 
Appendix A, the technique allows the decoupling of the SED 
of the underlying stellar population from that of the ionized 
gas. This step is needed to derive the intensities of the differ- 
ent emission lines detected in each spectrum accurately. In 
addition, when the signal-to-noise ratio is high enough this 
fitting technique can be used to derive the physical parame- 
ters that characterise the composite stellar-population: the 
luminosity-weighted age, metallicity and dust attenuation. 



5.2.1 Distribution of the stellar populations 

As shown in Appendix A, the results from simulations indi- 
cate that an individual spectrum does not have, in general, 
the required signal-to-noise ratio to derive accurate values of 



the physical parameters that characterise the composite stel- 
lar population. However, the average spectra in concentric 
annuli do possess the required signal-to-noise, up to ~ 120 
arcsec (~ 0.4p25 or ~ 5.4 Kpc) in linear-projected galacto- 
centric radii. 

In order to study the distribution of the properties 
of the stellar population we obtained azimuthally aver- 
aged spectra by co-adding all the spectra in the mosaic 
within successive rings of 5 arcsec, starting from the cen- 
tral emission peak of the galaxy. These azimuthally aver- 
aged spectra were then analysed using the fitting proce- 
dure described above. Fig. 8 shows the radial distributions 
of the luminosity-weighted ages (left-panel) and metallici- 
ties (right-panel), derived by this method. The adopted error 
bars correspond to the typical errors estimated from the sim- 
ulations (Appendix A) for a signal-to-noise ratio of ~ 100. 

It is important to stress here that the accuracy of the 
absolute scale of the ages is not as good as the relative one, 
since the former is much more model dependent than the 
second. The age distribution shows, in general, a gradient, 
with older stellar populations at the central regions of the 
galaxy and younger ones at the outer regions. However, the 
oldest stellar populations are not found in the centre of the 
galaxy, but in a circumnuclear ring at ~ 25 arcsec radius, 
with an inversion of the gradient at smaller radii. This region 
corresponds to the central kpc of the galaxy. At larger radii, 
the stellar population becomes steadily younger, at a rate of 
~ 0.20 Gyr/kpc, up to a radius of ~ 100 arcsec (~ 4.5 kpc), 
and then becomes stable with an age of ~ 1 Gyr, on average. 
On the other hand, the radial distribution of the metallicity 
of the stellar populations does not show such a clear trend, 
remaining nearly constant at any radii with a larger disper- 
sion in the values. Both distributions show a similar shape 
than the one presented by (MacArthur et al. 2009), using 
a similar technique over classical slit-spectroscopy data, al- 
tough their data show a larger dispersion. 

The luminosity-weighted age and metallicity, derived 
from the plotted values are ~ 8.2 Gyrs and Z ~ 0.01 
([Z/H]'^ —0.27), consistent with the values derived from the 
integrated spectrum Table 2. 

The gradient found in the stellar population of NGC 628 
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Figure 9. Left- panel: Reconstructed, continuum subtracted Ho narrow-band image of NGC 628 obtained after interpolating the aperture 
photometry extraction of the SINGS narrow-band image at the location of the mosaic fibres as in Fig. 4. Right-panel: Pure Ha emission 
lines map of NGC 628 calculated from the IFS mosaic in units of 10~^^ ergs^^ crn"^ arcsec^^. No correction for dust extinction was 
applied to the map. All morphological structures of the galaxy are completely reproduced in the IFS image. 



(with the older stellar populations in the center), is common 
in late-type galaxies. Gadotti & dos Anjos (2001) found that 
~ 60% of this kind of galaxies show a negative color gradi- 
ent in the U — B and B — V colors that it is most likely 
due to a gradient in the age of the stellar populations, with 
the younger populations in the outer regions. Differences 
between the colors in the central bulges (redder) and the 
outer disks (bluer) are also reported in other studies (e.g. 
Peletier & Balcells 1996). NGC 628 was included in the sam- 
ple studied by Gadotti & dos Anjos (2001), and they re- 
ported a negative color gradient in this galaxy up to a ra- 
dius of ~5 arcmin, qualitatively in agreement with the trend 
found here based on the analysis of our IFS dataset. 

As mentioned before, Ganda et al. (2007) presented a 
study based on the IFS observation of the central core of 
NGC 628, using SAURON. The total field of view of their 
observation was 33" x41" , with a spatial sampling of ~ 1 arc- 
sec. This FOV fits inside one of our single pointings, being 
^ 30 times smaller in area than the FOV of the full dataset 
presented in this article. Due to the dithering technique ap- 
plied at the central pointings, our final spatial resolution 
is just 2 times worse than their study. On the other hand, 
their spectral resolution is 2 times better than ours, with 
the penalty that their spectral coverage (4800 — 5300 A) is 
~ 7 times shorter than ours. Despite these differences, their 
results can be compared with the ones presented here for 
the inner regions. 

Fig. 6 of Ganda et al. (2007) shows, among other 
things, the spatial distribution of the age and metallicity of 
NGC 628 derived from their analysis. Despite the very small 
FOV of their observations and the more reduced wavelength 
range sampled by their data, the basic structure described 
above can be seen. In particular, the age distribution ap- 



pears flat in the inner 15 arcsec radius, with a ring of older 
stellar populations at ~ 20 arcsec. On the other hand, the 
structure in the metallicity shown by Ganda et al. (2007) 
is more clumpy. The azimuthally averaged radial profile of 
both quantities, as presented in Fig Ala of that article, 
shows a clearer picture, with a similar trend in both param- 
eters as shown in Fig. 8 for the same radii (although with 
an offset of the absolute values) . The decline in the age of 
the stellar population inferred from our data at larger radii 
is out of the FOV of the SAURON observation. 

Despite this qualitative agreement between both re- 
sults, there are significant differences in the quantitative 
parameters. The values of age and metallicity derived by 
Ganda et al. (2007) for the central regions are younger (~2 
Gyrs), and richer {Z ~0.0-0.2), than the ones found by our 
analysis. These differences are due to the different tech- 
nique used in both analysis. The analysis performed by 
Ganda et al. (2007) that we are comparing with, is based on 
a single stellar population fitting technique. Indeed, when we 
perform a similar analysis, adopting our fitting techinique to 
compare with SSP, instead of allowing to mix multiple ones, 
we derive an age of ~1.4 Gyrs, in a better agreement with 
the reported value. 

As noted by the authors, the SSP analysis does not re- 
flect perfectly well the complex starformation history of late- 
type galaxies. In particular, when a galaxy has undergone 
two separate bursts of star formation, the age derived by 
this analysis will be biased towards the one of the youngest 
stars, while the metallicity will be biased towards the on of 
the older population, as studied and described in detail by 
Serra & Trager (2007). Due to this effect, the authors refined 
their analysis of the stellar populations, constraining the 
metallicity by adopting the empirical relation between this 



16 S. F. Sanchez et al. 

Ha 




ri V, 



Aa (arcsec) 

Figure 10. Zoom of the central ±30" region of the pure H« 
emission lines map of NGC 628 calculated from the IFS shown in 
previous figure. The countour plots shows the continuum emission 
in the same region. The ionized ring described by Wakker &; Adler 
(1995), and confirmed by Ganda et al. (2006) and FathiOT, is 
clearly seen. 



value and the velocity dispersion, derived by Thomas et al. 
(2005). By adopting this constrain, they derive and age of 
7.943 Gyrs and a metallicity of [Z/H] =-0.498 for the stellar 
population of the NGC628 (Table 5 of Ganda et al. 2007), 
in a complete agreement with our results, hsted in Table 2. 



5.2.2 Distribution of the emission lines 

The ionized gas in any spiral galaxy exhibits a com- 
plex structure, morphologically associated with the star- 
forming regions located along the spiral arms. To study 
their properties, these regions have perviously been targeted 
by narrow-band (e.g. Kennicutt & Hodge 1980) and Fabry- 
Perot imaging (e.g., Fathi07). In many cases these narrow- 
band images catch more than one single line (e.g. Ha, the 
[Nil] AA6548,6584 doublet or the [Sll] AA6716,6731 dou- 
blet), reducing their potential for the study the basic pa- 
rameters of the ionized gas, or requiring certain assumptions 
about the actual line ratios included within the filter width 
(Kennicutt et al. 2008). IFS allows a much more refined de- 
coupling of the different lines, and the production of maps 
of each individual emission line. 

As in the case of the integrated spectra, in order to ex- 
tract any physical information from the data set, we need 
first to identify the detected emission lines of the ionized 
gas and to decouple their emission from the stellar contin- 
uum. We use for this purpose the clean IFS mosaic version 
of NGC 628 (introduced previously in Section 5.1), which 
is free of bad fibres, foreground objects, and includes only 
those regions with high-quality spectrophotometric calibra- 
tion, and fibres with an average flux along the whole spectral 
range greater than 10~^^ ergs~^ cm~^ A~^. 

After this (relatively conservative) flux threshold was 



Table 6. Number of fibres in different mosaic versions, the num- 
ber in the observed mosaic refers to the total number of spectra 
considering all the pointings of NGC 628. The clean mosaic num- 
ber corresponds to the remaining fibres after the flux threshold 
cut as explained in the text. The percentage stands for number 
of clean fibres with respect to the original mosaic. 

Observed mosaic Clean mosaic Percentage 



11094 



6949 



63% 



applied, the number of fibres remaining in the clean mosaic 
was somewhat reduced, as shown in Table 6. The final num- 
ber of spectra in the clean mosaic accounts for ~ 63% of 
the total number of fibres in the original mosaic. The SSP 
fitting method was then applied to each individual spectrum 
of the clean mosaic, obtaining a residual spectrum per each 
fibre, in which the intensities for each detected emission line 
were measured following the techniques described in detail 
in Section 5. The statistical uncertainty in the measurement 
of the line flux was calculated by propagating the error as- 
sociated to the multi-component fitting and considering the 
signal-to-noise of the spectral region. 

As a result of the process described above, a set of mea- 
sured emission line intensities (plus associated errors) was 
obtained for each observed spectrum of the final clean mo- 
saic of NGC 628. From these sets of emission line intensi- 
ties, emission line maps were created by interpolating the 
intensities derived for each individual line in each individ- 
ual spectrum, based on the position tables of the clean mo- 
saics, and correcting for the dithering overlapping effects 
when appropriated. The interpolation was performed using 
E3D, adopting a natural- neighbour, non- linear interpolation 
scheme, with a final scale of 1" /pixel in the resulting maps. 
Regions in the borders of the mosaics and/or large regions 
in between the mosaic without signal are prone to artifacts 
created by the interpolation scheme, special masks were cre- 
ated in order to deal with those regions. Further, as many 
of the derived maps are based on a reliable measurement 
of H/3, a fiux threshold mask was created for each mosaic 
for those regions with an integrated H/3 fiux per fibre be- 
low lO"^** ergs"^cm"^ (i.e. 0.2 x 10"^® ergs"^cm"^ per 
pixel), which corresponds to a detection limit of ~ 5a. This 
conservative limit has been adopted to avoid contaminating 
our results with low signal-to-noise data, and to grant the 
accuracy of the absolute spectrophotometry. 

Practically any emission line map (or a combination 
of them) can be constructed from the detected emission 
lines in the mosaic. In this section, we present a number 
of maps corresponding to the most important emission lines 
and derived quantities of interest ®. The right-panel of Fig. 9 
shows the pure Ha emission line intensity map in units 
of 10"^® erg s~^ cm~^ arcsec"'^, obtained by the proce- 
dure described previously. For a purely qualitative compar- 
ison purpose, the left-panel of the same figure shows the 
continuum-subtracted Ha narrow-band image of NGC 628 
extracted from the SINGS ancillary data by interpolating 
the aperture photometry extraction at the location of the 
mosaic fibres as described in Section 4.1. A visual inspec- 
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Figure 11. Dust attenuation map of the ionized gas {Av), de- 
rived by comparing the measured Hq/H/3 line ratio map with 
the expected value for case-B recombination, assuming a homo- 
geneous distribution of the electron density and temperature, and 
the attenuation law by Cardelli et al. (1989). The regions without 
accurate measurements of either emission linos were masked. 



tion of both maps shows that they match remarkably weU 
in terms of the morphology. As the dominant gas ionization 
mechanism giving rise to the emission line in this galaxy is 
due to hot (OB) stars (Section 5.1.2), the Ha line intensity 
map shown in Fig. 9 traces the star-forming regions in this 
galaxy. The map displays a distinctive spiral structure with 
multiple H II regions and H II complexes of different sizes and 
morphology along the spiral arms. The brightest sources are 
located at the outer regions of the galaxy, with a particular 
giant region located at (Aa, A5) ~ (-40,120) arcsec. From 
this map we can also note the presence of diffuse emission 
along the spiral arms and in the intra-arms regions. 

FathiOT has recently provided with a similar Ha map, 
with twice the field-of-view of our IFS data, on the basis of 
Fabry-Perot and complementary narrow-band images. The 
distribution shown by both maps are remarkable similar. 
They identified 376 HII regions, many of them within the 
field of view of our IFS data. A detail comparison of each 
of these individual sources, and an analysis of their spec- 
troscopic properties will be presented elsewhere (Rosales- 
Ortega in prep.). 

Their data confirm the presence of a central struc- 
ture in NGC 628, characterized by a ring of ionized gas at 
~15-20", already described by Wakker & Adler (1995) and 
Ganda et al. (2006) in their H/? and [OIII] emissin line dis- 
tributions. Figure 10 shows a zoom of the Ha distribution 
already shown in Figure 9 (right panel), together with a 
counterplot of the continuum emission in the V-band (fig. 4, 
right panel). The spatial resolution for this region is better 
than the one in the overall Mosaic, due to the dithering per- 
formed only in the central pointing. The slightly oval ring 
structure in the ionized gas is clearly identified in the Ha 
map, showing a remarkable similarity with the H/3 distribu- 
tion shown by Ganda et al. (2006). It is interesting to note 



here that despite the difference in the original sampling of 
the SAURON instrument (used by Ganda et al. 2006, 2007), 
and PPAK, the apparent resolution of both maps are very 
similar. This ring is detected, at different levels, in all of the 
emission line spices analyzed in the current article, following 
the same basic morphological pattern. 



5.2.3 Distribution of the dust attenuation 



Besides the intrinsic significance of the emission lines map 
presented previously, the most interesting application of the 
2D distribution of the different line intensities resides in the 
calculation of well-known line ratios that can be translated 
to physical properties of the emitting gas. With the aid of 
the IFS data we can investigate for the first time, the point- 
to-point variation of these physical properties over a consid- 
erable area on the surface of a galaxy. 

Fig. 11 shows the distribution of the dust extinction de- 
rived from the Ha/H/3 line ratio, in terms of the visual ex- 
tinction Av (in magnitudes) , assuming, as in the case of the 
analysis of the integrated spectrum, case-B recombination 
with an electron temperature of Te ~ 10* K, and adopt- 
ing the Cardelli et al. (1989) extinction law with Rv, the 
total to selective extinction ratio, equal to 3.1. Although 
other Balmer lines were detected at different locations in 
the galaxy. Ha and H/3 are obviously the lines with the 
highest signal-to-noise and therefore are the most adequate 
for the determination of Ay. The dust extinction map was 
derived for those locations where the intensity of the H/3 
line was above the adopted flux threshold. These regions 
clearly follow the spiral arm structure. The dust shows a 
clumpy distribution, a rich structure and large variations 
even within the same HII region or complex. There is no 
apparent trend of the extinction along the spiral arms or 
in any radial direction. The average extinction derived from 
the values shown in the map is Ay ^ 1.24± 0.76 mag, which 
is slightly larger but comparable to the value derived using 
the integrated spectrum {Ay ~ 1.04). The reason for this 
discrepancy may reside in the fact that the ionized gas com- 
ponent of the integrated spectrum is dominated by the spec- 
tra of the outer regions, where the intensity of the emission 
lines is stronger and the extinction is somewhat lower. Ex- 
tinctions greater than 2.5 are found in specific regions, e.g. 
the central zone of the giant HII complex at ~ (-40,120), 
while others are found in compact HII regions along both 
spiral arms. Previous studies have estimated the extinc- 
tion at different locations within this galaxy targeting in- 
dividual HII regions (e.g. McCall et al. 1985; Belley & Roy 
1992; Petersen & Gammelgaard 1996), the derived extinc- 
tions from these studies are consistent with the range of Ay 
values found in this work. 



5.2.4 Ionization conditions 

As discussed previously in Section 5.1.2, the source and 
structure of the ionization can be in principle investigated 
by exploring the line ratio maps of typical diagnostic in- 
dices. Fig. 12 shows on the left-panel the line ratio map 
of [Nil] A6584/Ha (N2Hc) and on the right-panel, the 
[Qui] A5007/H/3 (OsHp) ratio map, in logarithmic scale. 
Note that, given the proximity of the emission lines, these 
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Figure 12. 2D spatial properties of NGC628, classical diagnostic line ratio maps. Left-panel: [Nil] A6584/Ha line ratio map. Right 
panel: [OIII] A5007/H/3 line ratio map; both in logarithmic scale. 



indices are almost reddening-independent. The A^2-H'q shows 
a very homogeneous behaviour, with small variations in indi- 
vidual regions, some of them showing lower values in the cen- 
tral part of the regions (e.g the knot at Aa, A5 ~ 40,-130). 
The average ratio derived from this map is -0.55 ± 0.12 dex, 
which is excellent agreement with the value derived from 
the integrated spectrum {N2Ha = -0.56). However, some 
regions at the outer part of the galaxy show lower ratios 
consistent with values ~ —0.75, e.g {Aa,A5) ~ (-50,-160), 
(0,200), which might suggest a gradient of this index de- 
creasing from the inner to the outer regions of the galaxy. 
On the other hand, the O^Hp shows a clear gradient along 
the spiral arms with lower ratios towards the inner regions 
and greater values at the outer part of the galaxy. The aver- 
age value of this index derived from the map is -0.50 ± 0.25 
dex, which again is in good agreement with the integrated 
value of -0.48. It is interesting to note that the regions with 
the highest values of OsHp are coincident with the zones of 
the lowest N2Ha ratios. Given that the [Nil] emission origi- 
nates in the singly ionized regions, between the fully ionized 
and the partially ionized zones, the N2Ha ratio traces the 
changes in the local ionization, while the [O III] originates 
in the fully ionized zones, tracing the strength of the ioniza- 
tion. Therefore, the distribution found in these diagnostic 
maps may indicate that the ionization is stronger in the 
outer parts of the spiral arms, than in the central regions. 
The values of both ratios at any location in the galaxy are 
consistent with ionization produced by hot OB stars, as ex- 
pected. In particular, there is no evidence of ionization due 
to shocks and/or nuclear activity. 

An additional way of studying the 2D distribution of the 
galaxy properties consist in obtaining azimuthally-averaged 
radial spectra, from which radial average properties can be 
derived. Taking again as a base the clean IFS mosaic ver- 
sion of the galaxy, radial average spectra were obtained by 
co-adding all the spectra within successive rings of 10 arc- 
sec, starting from the central reference point. An average 



spectrum was obtained for each single annulus at a given 
radius. The radial average spectra were then analysed us- 
ing the same fitting procedures described before. Although 
the derived spectra present more signal-to-noise than the 
single-fibre case, the measured emission lines were corrected 
by extinction using only the Hq:/H/3 ratio for consistency 
with the previous method. 

Fig. 13 shows the BPT diagnostic diagram 
(Baldwin et al. 1981) for the N2HC and O3H13 line 
ratios. The bluish-asterisk symbols (in the online version) 
correspond to the derived values of the emission line ratios 
for 300 randomly selected spectra distributed uniformly 
over the surface covered by the emission line ratio maps 
shown in Fig. 12. The reddish diamonds correspond to 
the radial average values introduced above, while the red 
solid-square indicates the average line ratio derived from 
the integrated spectrum of the galaxy. The colour-coding 
of both the individual fibres and the radial-average samples 
is related to the spatial position of a given fibre/annulus. 
Lighter tones correspond to the inner regions of the galaxy, 
while darker colours correspond to positions with increasing 
galactocentric radius (arrow direction). The dark-thick 
line corresponds to the theoretical boundaries dividing 
the starburst region from other types of ionization using 
the parametrization provided by Kewley et al. (2001). 
The dashed lines represent the ±0.1 dex variation. As 
expected, the ionization is dominated by hot stars (OB 
class) associated with star-forming regions, at any location 
in the galaxy. A clear trend can be noticed in Fig. 13 
for both samples, the spectra corresponding to the inner 
regions tend to have lower line ratios for both indices; for 
regions at the outer part of the galaxy, the ratios increase 
approaching the theoretical boundary. Part of the reason 
for this behaviour is that inner parts of the galaxy lack 
emission in [O II] and [O III] , while towards the outer parts, 
the emission from these species is prominent, increasing the 
two line ratios involved in this diagram 
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Figure 13. N2HC, vs. O^Hp diagnostic diagram for NGC 628. 
The bluish symbols correspond to the derived values of the emis- 
sion line ratios for 300 randomly selected spectra distributed 
uniformly over the surface covered by the emission line ratio 
maps shown in Fig. 12. The reddish diamonds correspond to the 
azimuthally-averaged radial values. The red solid-square indicates 
the ratio derived from the integrated spectrum of the galaxy. 
Lighter tones correspond to inner regions of the galaxy, darker 
colours to increasing galactocentric distance, outlined by the ar- 
row. In the radial case, the size of the symbols increases with 
radius. 
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Figure 14. Distribution of the ionization parameter 

\ogu, in logarithm scale, derived from the dust corrected 
[Oil] A3727/[OIII] A5007 line ratio. 
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Figure 15. Surface star-formation rate (SFR) in units of 10~ 
Mq yr~^, derived from the Ha line intensity map shown in Fig. 
10, after correcting by dust attenuation. This map traces the dis- 
tribution of the ionized gas, not the young stars, which form much 
more compact structures than the ones shown in the figure. 



The differences in the hne ratios at different locations 
in the galaxy may be also driven by the strength of the ion- 
ization field, being stronger in the regions with higher N^Hc 
line ratio, and lower O^H^ one. In order to investigate this 
point. Fig. 14 shows the distribution of the ionization pa- 
rameter in logarithmic scale, logu, calculated accordingly 
to Equation 1 (Diaz & Perez-Montero 2000), which is based 
on the dust-corrected [O II] /[O III] ratio. Both lines were cor- 
rected by extinction using the Ay map discussed previously. 
The logu map shows that the ionization is indeed stronger 
in the outer parts of the spiral arms, as the previous maps 
suggested. However, we can also distinguish a good degree 
of ionization structure in individual regions within the spi- 
ral arms, with higher values of logu corresponding to the 
geometrical centroids of the HII regions, as one might ex- 
pect in the scenario of a central star cluster or association 
embedded within an HII region. Given that the dominant 
source of the ionization is radiation from OB stars, this re- 
sult indicates that the star-formation rate is stronger in the 
outer regions of the galaxy. A striking feature of this map 
is the presence of regions with low values of the ionization 
parameter (blue colour) located mainly in an specific region, 
corresponding to the first quadrant, north-west part of the 
galaxy. These low log u regions are found mainly at the edges 
of giant H II complexes, but they are also found as individ- 
ual regions. The nature of these low ionization regions will 
be investigated thereinafter. The mean value of log u de- 
rived from this map is -3.50 ± 0.26 dex, which is 0.08 dex 
higher than the value derived from the integrated spectrum 
of NGC 628. 



5.2.5 Surface star-formation rate 

The surface SFR can be derived using the emission line dis- 
tribution of Ha or [O II] A3727. However, the latter line 
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Figure 16. 2D distribution of the R23 and [O III]/[N II] mctallicity indices for NGC 628. Both maps show a gradient of lower values in 
the inner regions to higher values to the outer part of the galaxy. 



suffers from more inaccuracies in both the spectrophoto- 
metric cahbration and flux determination, since it hes just 
at the blue edge of the spectral range covered by our data 
set. Therefore, only the Ha intensity map was used in this 
analysis. Following our previous prescription to derive the 
integrated SFR of the galaxy, the Ha intensity map was 
corrected for attenuation, adopting the dust distribution de- 
scribed in Section 5.2.3, and the same attenuation law as 
in Section 5.1. Then, the dust corrected flux densities were 
transformed to absolute luminosities, using the same dis- 
tance and cosmology described before. These luminosities 
were used to obtain the surface SFR, by adopting the clas- 
sical relations described by Kennicutt (1998). 

Fig. 15 shows the derived distribution for the surface 
SFR, in units of 10^^ Mq yr^^ arcsec"'^. As expected, the 
stronger star-formation is distributed in clumps along the 
spiral arms, following a similar morphology to the distribu- 
tion of the Ha emission (i.e. the HII regions). In addition, 
the direct transformation between extinction corrected Ha 
flux and SFR produce an apparent, low level, star-formation 
rate in the inter-arms regions, two orders of magnitude lower 
than that of the HII regions. However, it is unclear how 
closely is associated this derived SFR with the inter-arm 
Ha emission. Ionizing photons leaking from the HII regions 
and stochastics in the conversion between both parameters 
may become important at this intensity level. As already 
shown by the distribution of the N2Ha and O^Hp line ra- 
tios (Fig. 12), and the ionization parameter (Fig. 14), there 
is an increase of the star-formation rate from the inner to 
the outer regions of the spiral arms. 

The global SFR can be derived by co-adding the surface 
SFR over the area shown in the flgure, obtaining a value of ~ 
1.1 Mq yr^'^. By applying the same aperture and sampling 
correction applied to the integrated Ha flux in Section 5.1, 
the estimated SFR is ^ 2.4 Mq yr~ , which is practically 
identical to the estimate derived based on the flux of Ha in 
the integrated spectrum. 



5.2.6 Distribution of the gas-phase oxygen abundance 

One of the main purposes of the PINGS dataset is to study 
the metallicity content of the galaxy in a 2D context. As dis- 
cussed previously, the oxygen abundance is normally derived 
using a variety of diagnostic methods based on strong lines, 
froiu which the most common one is the 7?23 indicator. The 
dust-corrected emission line maps of [O II] and [O III] pro- 
vide the means to explore the behavior of this indicator. The 
left panel of Fig. 16 shows the R23 map of NGC 628. From 
a flrst visual inspection, it is evident that the map presents 
a clear gradient of lower values of 7?23 in the inner regions 
of the galaxy, to higher values to the outer parts, specially 
along the spiral arms. The average R23 value derived from 
the map is 2.41 ± 1.37, which is in good agreement with the 
value inferred from the integrated spectrum {R23 = 2.57). 
The lower values of i?23 in the central regions of the galaxy 
are expected, given the low emission of [O II] and (especially) 
[O III] , as inferred in figure Fig. 13. For well-deflned HII re- 
gions, the value of i?23 seems to be constant in the majority 
of the cases; however, there are regions in which the value 
of i?23 varies within the HII region or complex, showing 
some level of structure. Considering the 0/H values derived 
after different calibrators from the integrated spectrum of 
NGC 628 in Section 5.1.3 (12 + log(0/H) ~ 8.7), and previ- 
ous abundance determinations from individual HII regions 
in this galaxy (e.g. McCall et al. 1985; Ferguson et al. 1998), 
suggest that the average oxygen abundance of NGC 628 cor- 
responds to the high metallicity regime. In this scenario, the 
double-valued nature of the i?23 index can be broken, and 
the gradient of higher-to-lower values of _R23 from the inner 
to the outer parts of the galaxy shown in this map would 
correspond to the decrement of the R2Z ratio in the upper- 
branch of a O/H vs. R23 diagram, and therefore, to a true 
metallicity gradient of the galaxy, as previous studies based 
on individual H II regions have shown. 

Another popular metallicity-sensitive index is given by 
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the [OIII] A5007/[NII] A6584 ratio. Pettini & Pagel (2004) 
suggested the use of this ratio (in a modified version) 
as abundance indicator suitable for the analysis of high- 
redshift galaxies. However, the direct application of this in- 
dicator may be too simplistic, the reasons being that the 
[O III] /[Nil] ratio is strongly dependent on the excitation 
of the nebula, and it is also sensitive to both the ioniza- 
tion parameter and to the age of the cluster of exciting 
stars (Dopita et al. 2000). As suggested by Dopita et al. 
(2006), the [OIII] A5007/[NII] A6584 ratio can only be used 
as an abundance indicator only when a characteristic age 
of the exciting clusters can be assumed, and its applica- 
tion in single galaxies should be taken with caution. How- 
ever, the combination of this ratio with, for example, the 
[Nil] A6584/[0II] A3727 ratio, can provide a good diag- 
nostic of both metallicity and age of the ionizing source 
(Dopita et al. 2006). The right panel of Fig. 16 shows the 
map of the [OIII] A5007/[NII] A6584 ratio of NGC 628 in 
logarithmic scale. It is worthy to note that this ratio spans 
for more than two orders of magnitude. Similarly to the 
7?23 map, this panel shows a gradient of lower values of the 
[O III]/[N II] ratio from the inner part of the galaxy, to higher 
values in the outer parts. The average value derived from the 
map is -0.35 ± 0.36 dex, compared to -0.40 obtained from 
the integrated spectrum. The map shows a smoother distri- 
bution within individual regions compared to the _R23 map. 
With the corresponding cautions considering the known de- 
pendences of this index with the functional parameters of the 
HII regions, the [OIII] /[Nil] gradient found for NGC 628 
would correspond to a metallicity gradient of the galaxy, 
with some level of inhomogeneity in individual regions as 
shown in this map. In order to determine the existence of 
abundance variations within the same region of the galaxy, 
it would be required to co-add the spectra corresponding to 
regions with similar ionization conditions, and then perform 
an analysis based on single regions. This detailed analysis, 
together with a complete study of the radial abundance gra- 
dient of this galaxy will be presented in a forthcoming paper 
(Rosales-Ortega et al. in preparation). 



In Fig. 17 we present the oxygen abundance 2D map 
of NGC 628 based on the 03N2 calibrator in units of 
12 + log(0/H), obtained by the manipulation of the emis- 
sion line maps presented before. No other 2D abundance 
maps were attempted since the parametrizations of the dif- 
ferent calibrators do not allow a simple image manipulation 
(large blank areas produce convergence problems which lim- 
its the application of any interpolation routine). However, as 
discussed previously in Section 5.1.3, the 03N2 abundance 
determination follows qualitatively the same pattern as the 
other more elaborated calibrators based on R23 and the ff- 
Te relation. The map of Fig. 17 shows a clear gradient of 
higher oxygen metallicity values from the inner part to the 
outer part of the galaxy, and along the spiral arms. The 
most metal deficient regions are found at the extremes of 
the spiral arms, a prominent region is found at {Aa,AS) 
~ (-40,-60). The average value of the oxygen abundance 
derived from this map is 8.75 ±0.17, which is equivalent to 
the oxygen abundance derived from the integrated spectrum 
using the same calibrator, i.e. 8.71 ±0.11. 
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Figure 17. Oxygen abundance map of NGC 628 derived by ap- 
plying the 03N2 calibrator (Pettini & Pagel 2004) to the emis- 
sion line maps of the galaxy. The figure shows a clear gradient in 
metallicity, with more abundant regions in the inner part or the 
galaxy. 



Table 7. Summary of the main properties of the galaxy. 



Parameter 


Integrated 


Average 


Age (Gyr) 


8.95 




[Z/H] 


-0.44 




Av stellar (mag) 


0.4 




Av (mag) 


1.04 


1.24±0.76 


Iog([NII] A6584/Ha) 


-0.56 


-0.55 ±0.12 


log([OIII] A5007/H/3) 


-0.48 


-0.50 ±0.25 


logM 


-3.58 


-3.50 ±0.26 


SFR (M0 yr-l) 


2.4 


2.4* 


log -R23 


0.41 


0.38±0.13 


log([OIII] A5007/[NII] A6584) 


-0.40 


-0.35 ±0.36 


12 + log(0/H)o3N2 


8.71±0.11 


8.75±0.17 


12 + log(0/H)t 


8.69 ±0.31 





* Co-added. 

f Mean abundance among the various empirical calibrators. 

6 DISCUSSION AND CONCLUSIONS 

In this article we presented wide-field IFS observations of a 
substantial fraction of the nearby galaxy NGC 628 (~ 70% of 
its projected optical size). The observations were performed 
by adopting a mosaicking procedure, using the PPAK IFU 
(mounted at the 3.5m telescope at Calar Alto), which pro- 
vides a sampling of 2.7 arcsec and a FOV larger than 1 ar- 
cmin (diameter). The total dataset comprises 11094 individ- 
ual spectra, covering a nearly circular FOV of '-^ 6 arcmin in 
diameter, centered on the bulge of the galaxy. To our knowl- 
edge, this is the largest spectroscopic survey ever made of a 
single nearby galaxy. A detailed flux calibration technique 
has been applied, which granted a spectrophotometric accu- 
racy of ~ 0.2 mag, down to a surface brightness of ~ 22 mag 
arcsec"^, comprising 6949 spectra from the original dataset. 
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The spectrophotometrically calibrated data have been 
analysed both as a single integrated spectrum, which char- 
acterises the global properties of the galaxy, and using 
each individual spectrum to determine the spatial variation 
of galaxy properties. The stellar and ionized gas compo- 
nents have been decoupled by a fitting technique which uses 
SSP templates and single Gaussian functions respectively 
to characterise these components. The main properties ex- 
tracted from the analysis of the stellar continuum emission 
were the luminosity-weighted ages, metallicities and dust at- 
tenuation of the composite stellar population. On the other 
hand, for the gas content, we derived the morphology of the 
ionized gas, the dust content, the spatially resolved SFR, 
the source of the ionization, and the oxygen abundance. 

Table 7 lists the main properties derived in this study 
for NGC 628. For each of the estimated parameters we list 
the values derived based on the analysis of the integrated 
spectrum (labelled as Integrated), and the average values 
derived from the 2D maps (labelled as Average). When 
more than one technique was used to derive the values, we 
listed the values derived using the same technique for both 
datasets. In all cases, the derived values are compatible. 

Regarding the spatial distribution of the different anal- 
ysed properties, a gradient was found in the luminosity- 
weighted age of the stellar populations from the central re- 
gions, having older stellar populations, to the outer regions, 
having younger components. On the other hand, the radial 
dependence of the stellar metallicity is weaker, showing a 
slight or no gradient. If there is any gradient, the richer 
stellar populations are located at the central regions, while 
the poorer ones are at the outer regions. All together these 
results are consistent with our current understanding of 
how the late-type galaxies grow along cosmological epochs, 
that is inside-out (Barden et al. 2005). The older and more 
metal-rich stellar populations, formed before, on the basis 
of recycled gas from previous populations, are located in 
the center, while the younger and less metal rich stars are 
in the outer regions. It is interesting to note that the oxy- 
gen abundance in the ionized gas follows nearly the same 
distribution (Fig. 17), with larger abundances in the inner 
regions and lower abundances in the outer ones. This gra- 
dient was previously reported by McCall et al. (1985), and 
confirmed by Belley & Roy (1992), van Zee et al. (1998) and 
Ferguson et al. (1998). 

In the case of the stellar populations, there seems to 
be an inversion of the age gradient at the very centre of 
the galaxy (~ 1 kpc) , where apparently there exist a ring of 
old stars at this distance, with a trend to younger ones at 
the very center. These results were previously reported by 
Ganda et al. (2007), on the basis of their IFS observations 
at the very centre of this galaxy. Similar results have been 
reported in other galaxies, mostly Sa/SO, where the inner 
regions of their bulges present bluer colors, consistent with 
younger stellar populations (e.g. Deharveng et al. 1997). Ac- 
cording to this scenario, the circumnuclear star-formation 
rings detected in Sa-Sbc galaxies might show evidence of 
gas transfer due to radial motions in these galaxies (e.g. 
Knapen et al. 2006). 

Spatially-resolved maps of the emission line intensities 
and physical properties were derived for NGC 628. Contrary 
to previous attempts to perform a 2D wide-field analysis 
based on narrow-band (or Fabry-Perot) imaging, which only 



allowed a basic analysis of the physical parameters and/or 
required assumptions on the line ratios included within in- 
dividual filters (e.g. Ha), the emission line maps presented 
in this paper were constructed from individual (deblended) 
emission fines at any discrete spatial location of the galaxy, 
where enough signal-to-noise was found. This fact allowed 
to investigate the point-to-point variation of the physical 
properties over a considerable area on the galaxy. Extinc- 
tion, ionization, and metallicity-sensitive indicators maps 
were derived from reddening corrected emission line maps. 
In general, they show that the ionized gas in these spiral 
galaxies exhibits a complex structure, morphologically asso- 
ciated with the star forming regions located along the spiral 
arms. 

In general, it is found that the dominant ionization pro- 
cess is due to the OB hot-stars associated with star-forming 
regions, based on the distribution of the classical diagnos- 
tic lines. The ionization is stronger along the spiral arms, 
associated with the HII regions, and more intense in the 
outer than in the inner ones. Indeed, the surface SFR is 
an order of magnitude stronger in the outer HII regions, 
at distance larger than ~ 100 arcsec (4.5 kpc), than in the 
inner ones. Considering that in these outer regions there 
is a lower mass density, the growing rate of stellar mass is 
considerably larger there than in the inner ones. Therefore, 
even taking into account the circumnuclear star-formation 
rate, the growth of the galaxy is dominated by the inside-out 
process. 

There is growing evidence of a central structure in NGC 
628, widely discussed by Fathi07. This structure is char- 
acterized by a ring of ionized gas (Wakker & Adler 1995; 
Ganda et al. 2006; Fathi et al. 2007), and a drop in the 
B-V colors (Natafi et al. 1992; Fathi et al. 2007). The kine- 
matics analysis of the stellar and gas component within this 
ring, in comparison with that of the rest of the galaxy, in- 
dicates that there is a rapidly rotating component, most 
probably a dynamically cold disc, which could be built by 
inflow from the outer parts towards the center (Daigle et al. 
2006; Ganda et al. 2006; Fathi et al. 2007). Fathi07 found 
that this circumnuclear ring is that the expected location 
of the inner Lindblad resonance radius, where the gas is in- 
deed expected to accumulate, due to non-circular motions 
exerted by a bar or spiral arms. This gas concentration may 
lead to a low luminosity starburst (e.g., Knapen et al. 1995). 
Our results show, for the first time, that (1) the gas in the 
circumnuclear ring is ionized by starformation, and (2) the 
ring is spatially coincident with a decrease in the luminosity- 
weighted ages of the stellar populations, which causes the 
previously observed drop in the B — V colors. All together, 
they support the interpretation outlined by Fathi07. 

The integrated spectrum of the galaxy was used to anal- 
yse the chemical abundance of the galaxy by means of dif- 
ferent empirical diagnostic abundance calibrators. We con- 
firmed the validity of the determination of the integrated 
chemical composition of a galaxy through the analysis of 
the global-emission line spectra, as previously found by 
other authors (Kobulnicky et al. 1999; Pilyugin et al. 2004; 
Moustakas & Kennicutt 2006a), under the assumption that 
the observed emission lines arise via photoionization from 
young, massive stars. From the set of calibrators used for 
this purpose, the R23 methods based on photoionization 
modelling provide higher values of the oxygen abundance 
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(i.e. M91 & KK04), followed by those methods which con- 
sider one of two emission-line indices (N2, 03N2) and the 
ff-Te method. The differences between different calibrators 
are as large as ~ ±0.7 dex, especially for the high values of 
0/H derived from the R23 methods compared to the other 
ones. If an arbitrary negative offset is applied (as suggested 
by the literature), the 0/H values seem to be consistent 
for an average integrated oxygen abundance of NGC 628 of 
12 4- log(0/H) ~ 8.7, in agreement with previous studies. 
Note that the validity of the abundance determination de- 
pends mainly on the chosen calibrator used to derive the 
chemical composition, and to a second order, on the SSP 
fitting decoupling. 

The spatially resolved distribution of the abundance 
shows a clear gradient of higher oxygen metallicity values 
from the inner part to the outer part of the galaxy, and 
along the spiral arms Fig. 17. However, in some instances, 
the value of the oxygen abundance (and other physical prop- 
erties like extinction and the ionization parameter) varies 
within what would be considered a classical well-defined 
HII region (or HII complex), showing some level of struc- 
ture. Indeed, the 2D character of the PINGS data allow us 
to study the small-scale variation of the spectra within a 
given emitting area. As suggested in Paper I, the values of 
the emission line ratios measured using different extraction 
apertures vary considerably as a function of the aperture 
size, and that the scatter of the central value is larger than 
the statistical error in the measurements, reflecting that this 
might in fact be a physical effect. If we take as a premise 
that for a sufficiently large HII region, the emission line 
measurements are aperture and spatial dependent, i.e. that 
the light is emitted under different physical conditions, by 
gas in different degrees of ionization, and modified by differ- 
ent amounts of reddening (and therefore providing different 
elemental ionic abundances), the level of structure seen in 
the 2D maps of NGC 628 may be due to the intrinsic distri- 
bution of the ionizing sources, gas content, dust extinction 
and ionization structure within a given region, i.e. we would 
be sampling real point-to-point variations of the physical 
properties within a H II region. 

The emission line maps presented in this paper proved 
to be useful in describing the general 2D properties of the 
selected galaxies. However, the conclusions raised from them 
were based on general trends that depend, to a certain level, 
on the interpolation scheme applied in order to derive the 
pixel-resolved maps. More robust conclusions can only be 
drawn by analysing specific individual regions within the 
given galaxy, or by co-adding spectra of regions with the 
same physical properties and comparing the results in the 
2D context. This approach will be followed in subsequent 
papers in order to study the spatially-resolved spectroscopic 
properties of this and the rest of the PINGS galaxies. 
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APPENDIX A: MULTIPLE STELLAR 
POPULATION FITTING 

A substantial fraction of the stellar continuum emission in 
early-type galaxies is dominated by a single-stellar popula- 
tion, reflecting the fact that their cosmological evolution is, 
in many cases, well described by a monolithic collapse model. 
In this scenario most of the star formation took place in 
these galaxies for a short time at early times, and therefore, 
their SFH can be well described by a single star-formation 
process. As a consequence, their stellar continuum can be 
easily described by a single-stellar population. 

It is well known that the spectral energy distribution 
(SED) of simple stellar populations (chemically homoge- 
neous and coeval stellar systems) depends on a set of first 
principles (e.g. initial mass function, star formation rate, 
stellar isochrones, metallicity, etc), from which it is possi- 
ble to generate the spectra of synthetic stellar populations. 
This technique, known as evolutionary synthesis modeling 
(e.g. Tinsley 1980), has been widely used to unveil the stellar 
population content of galaxies by reconciling the observed 
spectral energy distributions with those predicted by the 
theoretical framework. Unfortunately the variation of dif- 
ferent physical quantities governing the evolution of stellar 
populations produce similar effects in the integrated light 
of those systems, leading to a situation in which the ob- 
servational data is affected by undesirable degeneracies, like 
the widely mentioned one between age and metallicity (e.g. 
Oconnell 1976; Aaronson et al. 1978; Worthey 1994). How- 
ever, the use of spectrophotometrically calibrated spectra 
and the sampling of a wide wavelength range, as in the 
present data set, helps to break the degeneracy and al- 
lows the derivation reliable physical parameters by fitting 
the full spectral distribution with single-stellar populations 
(Cardiel et al. 2003). 

However, the simple assumption that a single-stellar 
population describes well the SED of a galaxy is not valid 
for late-type galaxies. These objects present complex star 
formation histories, with different episodes of activity, of 
variable intensity and time scale. Therefore, a single-stellar 
population does not reproduce well their stellar emission, 
in general terms. A different technique, known as full- 
spectrum modeling, involving the linear combination of mul- 
tiple stellar populations and the non-linear effects of dust 
attenuation, has been developed to reconstruct their stel- 
lar populations (e.g. Fernandes et al. 2005; Ocvirk et al. 
2006; Sarzi et al. 2006; Koleva et al. 2009; MacArthur et al. 
2009). 

In general, these reconstructions require a wide wave- 
length range to probe the hot, young stars and the cool. 
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old stars simultaneously. They also require the best spec- 
trophotometric calibration to disentangle the effects of age, 
metallicity and dust attenuation. Although different imple- 
mentations of this technique have some differences, they 
are very similar in their basis, as described above. The ex- 
tracted information from the multi stellar-population mod- 
elling differs in the different implementations. In some cases 
the luminosity (or mass) weight ages and metallicities are 
derived, based on the linear combination of different models 
(e.g. Sarzi et al. 2006). In other implementations the frac- 
tion of light (or mass) of different stellar populations (e.g. 
Stoklasova et al. 2009), or the fraction of light (or mass) 
that corresponds to old or young stellar populations (e.g. 
Mac Arthur et al. 2009) are derived. 

For our particular analysis, it is necessary to model ac- 
curately the underlying continuum not only to understand 
the nature of the stellar population, but also to obtain a 
good representation of the continuum in order to decou- 
ple it from the emission lines produced by the ionized gas. 
Therefore, even in the case that the combination of SSPs is 
strongly degenerate, and the created model has no physical 
meaning, it could be partially useful for our purposes. 

The implementation of the multi stellar-population fit- 
ting technique used in this article is part of the FIT3D 
package (e.g. Sanchez et al. 2007b, c), and it has been re- 
cently used by (e.g. Stoklasova et al. 2009) in the analysis 
of the stellar population of nearby galaxies based on IPS 
data. This package enables linear fits of a combination of 
SSPs, and non-linear ones of emission-lines plus an underly- 
ing stellar population. It also includes routines to extract the 
2D distribution of the different derived parameters based on 
the analysis of the stellar population (age, metallicity, dust) 
and/or the analysis of the emission lines (flux intensity, sys- 
temic velocity and velocity dispersion). The package will be 
delivered to the public in a forthcoming paper. The basic 
steps of the fitting algorithm, spectrum by spectrum, are 
the following: 



N -M 



X^^i, 



where 



Ri = wA Gi — 2_^ o-jF, 






J=l 



In the above expression, Wi is the weight of the i"* pixel, 
defined as 



and aj is the coefficient of the j* template in the final mod- 
elled spectrum, and M is the number of templates consid- 
ered in the fitting procedure. 

(v) Determine for which templates the fitting procedure 
produces negative coefficients in the linear combination (i.e. 
aj < 0). These templates will be excluded from the next 
iteration of the fitting procedure, that will be resumed in the 
step (iii). At each iteration, Af is decreased by the amount of 
excluded templates. This loop ends once all the coefficients 
are positive. 

(vi) The modified x^ corresponding to the considered 
dust attenuation is stored as Xay ' ^'^'^ ^^^ fitting procedure 
is repeated again starting at the step (iii), modifying the 
considered amount of dust and starting the procedure with 
the full set of templates (i.e. M = Mq, again). This loop ends 
once the dust attenuation has covered a pre-defined range 
of possible values, that it is defined as an input parameter. 

(vii) The best linear combination of templates and dust 
attenuations are selected based on the minimization of the 
stored x^ parameters. The final modelled spectrum fiux at 
Xi is given by 



(i) Read the input spectrum, determine the areas to be 
masked, and construct the variance map. Lets define d as 
the observed galaxy flux at wavelength Xi, and A'^ the num- 
ber of elements of the masked spectrum. 

(ii) Read the set of SSP template spectra, shift them to 
the systemic velocity of the considered spectrum, convolve 
them to match its spectral resolution and velocity disper- 
sion, and resample them to its wavelength solution. Lets 
define Fji as the flux of the i*'' wavelength of the j*'' tem- 
plate (once shifted, convolved and resampled), where Mo is 
the total number of considered templates. 

(iii) Apply a certain dust attenuation to the templates. 
The attenuation law of Cardelli, Clayton & Mathis (1989) 
was adopted, with a ratio of total to selective attenuation 
of Rv =3.1 (Jenkins 1987). Lets deflne F;^^ as the flux of 
the i' wavelength of the j' template, after applying the 
dust attenuation corresponding to a certain attenuation of 
Av magnitudes. 

(iv) Perform a linear least-square fltting of the input spec- 
trum with the set of SSP templates, using a modified x^ ^s 
a merit function to be minimised, with the form: 



S. = Xa,K 



j=i 



where in this case, only the M templates with positive coef- 
ficients aj are considered, and Av corresponds to that value 
which minimises Xa^ ■ 

(viii) The luminosity-weighted age (H) and metallicity 
(Z) of the underlying stellar population is then derived by 
the formulae: 



i/ = 5^a,//, 



and 



where Hj and Zj are the corresponding age and metallicity 
of the j^^ SSP template. These luminosity-weighted ages 
and metallicities should be considered as the luminosity- 
weighted a,ges and metallicities of the modeled stellar popula- 
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Figure Al. Spectra of the 72 SSP templates used in the population synthesis fits. Metallicity increases from the top-left to the 
bottom-right. Different ages are represented by different colours. All the spectra are normalized to the flux at 5000 A. 
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Figure A2. Results from the simulations. Each panel shows, from left to right, the comparison between the input luminosity-weighted 
age, metallicity and dust attenuation, with the recovered values by our fitting technique. The results are presented for different signal- 
to-noise ratios, decreasing from the upper to the lower panels. A random value of ±0.08 mag has been added to output attenuation to 
show all the values, since they are originally quantized to steps of 0.1 mag. 



tion, since they would match with the corresponding values 
if the population was composed by a single SSP. 

This algorithm shares a basis with many other proce- 
dures described in the literature (e.g. Fernandes et al. 2005, 
Sarzi et al. 2006, MacArthur et al. 2008). If required, low- 
frequencies in the spectrum can be fitted by adding to the 
fitted templates a polynomial function, or by multiplying 
them by that function. There are different origins for these 



low-frequencies in the spectrum, mostly related with defects 
in the reduction process. In the particular case of the data 
analyzed in the current article, the inclusion of this addi- 
tional step in the fitting procedure does not improve consid- 
erably the results. 

By construction, the fitting algorithm is useful for mask- 
ing undesirable regions in the spectrum. In general, it is 
needed to mask: (i) strong and variable night sky-line re- 
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gions, (ii) regions affected by defects in the CCD (like dead 
columns) whose effect was not completely removed during 
the data reduction process, (iii) regions affected by telluric 
absorption, not completely corrected during the flux calibra- 
tion process, and (iv) regions containing emission lines from 
ionized gas. 

As indicated above, one of the main goals of this fitting 
procedure is to provide with and accurate modelling of the 
underlying stellar population, and to derive a pure-emission 
line spectrum, given by 



Al Modelling the emission lines 

To derive the intensity of each detected emission line, the 
clean spectrum (d) is fitted to a single Gaussian func- 
tion per emission line plus a low order polynomial function. 
This fit is performed using non-linear fitting routines imple- 
mented in FIT3D, that were described in previous articles. 
Instead of fitting all the wavelength range in a row, we ex- 
tracted for each spectrum shorter wavelength ranges that 
sampled one or a few of the analysed emission lines, in or- 
der to characterise the residual continuum with the most 
simple polynomial function, and to simplify the fitting pro- 
cedure. When more than one emission was fitted simulta- 
neously, their systemic velocities and FWHMs were forced 
to be equal (since the FWHM is dominated by the spectral 



resolution) , in order to decrease the number of free param- 
eters and increase the accuracy of the deblending process 
(when required). These procedures are frequently used, and 
therefore, we will not explain them in more detail. 

A modelled gas-emission spectrum is created, based on 
the results of the last fitting procedure, using only the com- 
bination of Gaussian functions. This spectrum is given by 



E, = Y.B>^*GAUSS,ki\k,ak), 



where L is the number of emission lines considered in the 
model, Bfe is the integrated fiux of the fc*'' emission line, 
and GAUSSik is the corresponding normalized Gaussian 
function (with central wavelength Afe, and dispersion at), 
evaluated at the i*'' pixel. 

Finally, the original spectrum (Gi) is decontaminated 
by the gas emission, subtracting this modeled gas-emission 
spectrum, deriving a gas-free spectrum given by 



GF, = G^-E,. 

This spectrum is then used to model the stellar popula- 
tion, by applying the procedure described before, but with- 
out masking the emission lines spectral regions. We adopted 
the results from this second iteration as a better modeling 
of the stellar population. 
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A2 SSP template library 

As noted by MacArthur et al. (2009), this Idnd of analysis 
is always limited by the template library, which comprises a 
discrete sampling of the SSP ages and metallicities. It would 
be desired for the stellar library to be as complete as pos- 
sible, and non-redundant. However, this would require an 
exact match between the models and the data, which is not 
possible to achieve in general terms, in particular if the stel- 
lar population comprises more than one SSP. As a suitable 
solution for the analysis presented in this article, we adopted 
a SSP template library that covers the widest possible range 
of ages and metallicities. This template has already been 
demonstrated as useful for modelling the stellar population 
of a varied set of galaxies (see Sanchez et al. 2007c). 

The SSP models were created using the GISSEL 
code (Bruzual & Chariot 2003), assuming a Salpeter IMF 
(Salpeter 1955) (the use of Chabrier 2003, will not modify 
the results), for different ages and metallicities. We create 
72 models covering a discrete grid of 12 ages (5 Myr, 25 
Myr, 100 Myr, 290 Myr, 640 Myr, 0.9 Gyr, 1.4 Gyr, 2.5 
Gyr, 5 Gyr, 11 Gyr, 13 Gyr, and 17 Gyr), and 6 metallic- 
ities {Z =0.0001, 0.0004, 0.004, 0.008, 0.02 and 0.05). Fig- 
ure Al shows the spectra of all the templates, grouped by 
their metallicity, and normalized to their flux at 5000 A. 
The adopted library is similar in many respects to the one 
proposed by MacArthur et al. (2009). Both libraries have a 
similar number of templates, although the latter library has 
a wider coverage of the ages and a more reduced coverage 
of the metallicities of the stellar populations. 

There is a number of caveats when applying model SSPs 
to the integrated light of a star forming galaxy, which have 
been clearly identified by MacArthur et al. (2009). The most 
important one is to assume that the parameter space covered 
by the empirical library represents well that of the real data. 
However in general, libraries are based on stars in the solar 
neighborhood, and therefore it is not granted that they rep- 
resent well the stellar populations in other galaxies (or even 
in other regions of our Galaxy). There are other potential 
problems related to the particular selected templates, since 
it is well known that the Bruzual& Chariot (2003) models 
have problems when dealing with the non-solar abundance 
ratios. Most of these problems are not particularly impor- 
tant in the context of our science case, since (i) our primary 
goal is to model the stellar population to analyse the ionized 
gas emission and (ii) we are interested in relative assessments 
about the nature of the stellar populations in different re- 
gions of the studied galaxy, and not in absolute values. 

In addition, it is important to note that the treatment 
of the dust attenuation may affect the resulting derived pa- 
rameters (i.e. the luminosity-weighted age and metallicity of 
the stellar population). In this particular implementation of 
the analysis we adopted the Cardelli et al. (1989) attenua- 
tion law, which may not be the optimal solution to study 
the dust attenuation in star-forming galaxies (e.g. Calzetti 
2001). MacArthur et al. (2009) adopted a completely dif- 
ferent attenuation law, based on the two-components dust 
model of Chariot & Fall (2000), which is particularly devel- 
oped to model the dust attenuation in star forming galaxies. 
Despite the conceptual differences between the two attenu- 
ation laws, their actual shapes are very similar in the wave- 
length range covered by our data. The final range of dust 



Table Al. Results of the simulations. 



S/N 


rms 


AAge/Age 


AZ 


AAv 


734.60 


0.14 


0.02 ±0.14 


0.0001 ±0.0014 


0.00 ±0.05 


90.28 


0.11 


-0.03 ±0.24 


0.0006 ±0.0047 


0.01 ±0.08 


43.74 


0.11 


-0.08 ±0.30 


0.0004 ±0.0074 


0.03 ±0.13 


8.60 


0.12 


-0.18 ±0.46 


0.0014 ±0.0126 


0.04 ±0.38 


1.69 


0.12 


-0.18 ±0.59 


0.0006 ±0.0164 


-0.06 ±0.51 


*82.77 


0.12 


0.21 ±0.27 


0.0010 ±0.0070 


-0.01 ±0.10 


"81.97 


0.13 


0.12 ±0.47 


0.0023 ±0.0075 


-0.01 ±0.23 



* Results when fitting the simulated spectra with a different li- 
brary of SSPs than the one used to create the simulation. 
** Results when the input simulated spectra consists of a single 
SSP, instead of a combination of them. 



attenuations explored by our fitting technique comprises Ay 
values up to 1 mag, with steps of 0.1 mag. 



A3 Accuracy of the derived parameters 

As mentioned before, the basic parameters derived by the 
analysis of the stellar component are the luminosity-weighted 
age, metallicity and dust attenuation of the composite stellar 
population. In order to assess the accuracy of these parame- 
ters we have performed a set of simulations. The simulations 
were performed while trying to match as closely as possi- 
ble the original data, especially the noise pattern, which is 
composed of both white noise corresponding to the photon- 
noise of the source and the background, and electronic noise 
from the detector; and non- white noise corresponding to 
defects/inaccuracies in the sky-subtraction, uncorrected de- 
fects in the CCD, etc. These noise patterns are different 
spectrum-to-spectrum, and wavelength-to-wavelength, and 
are clearly difficult to simulate on a simple analytical basis. 

Five different sets of simulations were performed. In 
each one, 500 simulated spectra were created with the same 
normalized flux at 5000 A, corresponding to 100, 10, 5, 1 
and 0.2 10~^® erg s~^ cm~^ A~^, respectively. For each sim- 
ulation, a composite stellar population was considered, com- 
prising 4 different SSP extracted randomly from the consid- 
ered library template, with a different relative contribution 
each one to the total flux: 53%, 27%, 13% and 6% at 5000 
A, respectively. This procedure was adopted to resemble in 
a simple way a possible star-forming history, with a domi- 
nant stellar population, and three different bursts or merg- 
ing events. Different experiments adopting a larger number 
of intermix stellar populations, and relative contributions of 
each one to the total were considered, with no significant 
variations on the derived conclusions. 

Once a particular spectrum was created, it was con- 
volved with the instrumental resolution of the observed data, 
and the corresponding noise pattern was added. The noise 
pattern was created on the basis of the residuals of the fit- 
ting procedure applied over the observational data. For each 
spectrum of the original mosaic, the fitting procedure cre- 
ates a residual spectrum, once the derived composite stellar 
population and emission line models were subtracted. This 
noise spectrum includes all the effects of the white noise (for 
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Figure A3. Results from the simulations. Each panel shows, from left to right, the comparison between the input luminosity-weighted 
age, metallicity and dust attenuation, with the recovered values by our fitting technique. The difference from Fig. A2 is that in this 
particular case, a different library of SSPs was adopted to fit the spectra, compared to the one used to simulate them, as described in 
the text. 
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Figure A4. Results from the simulations. Each panel shows, from left to right, the comparison between the input luminosity-weighted 
age, metallicity and dust attenuation, with the recovered values by our fitting technique. The difference from Fig. A2 is that in this 
particular case, the input model consists of a single SSPs, although it is fitted following the prescriptions shown in the text. 



a particular flux level), and non- white noise (deficient sky 
subtraction, CCD defects not corrected by the reduction, 
imperfections of the modelling process, etc.). From this set 
of noise spectra, one was randomly selected for each sim- 
ulated spectrum, fulfilUng the requirement that the simu- 
lated and original fluxes at ~ 5000 A match within a range 
of ~ 20%. The selected noise spectrum was then added to 
the simulated spectrum. 

The fitting procedure was then applied to the simulated 
spectra using the same steps adopted to analyse the observa- 
tional data, deriving, for each one, the luminosity-weighted 
ages, metallicities and dust attenuation, and a new residual 
spectrum. Table Al lists, for each set of simulated spectra, 
their average signal-to-noise ratio, the rms of the relative dif- 
ference between the input spectra (without noise) and the 
model derived by the fitting procedure and the difference 
between the input and the recovered parameters that define 
these models (luminosity-weighted age, metallicity and dust 
attenuation), and their corresponding rms. As expected the 
parameters are better recovered for the spectra with higher 
signal-to-noise ratios, requiring at least a S/N of ~ 50 per 
pixel to derive accurate results. Figure A2 illustrates these 




50 100 

Radius (arcsec) 

Figure A5. Radial distribution of the average S/N at ~ 5000 A 
of the spectra within the flux calibrated mosaic of NGC 628. 

results, showing for each set of simulations the comparison 
between the input and recovered parameters for each indi- 
vidual spectrum in the dataset. 
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A possible source of error in the interpretation of these 
simulations is the fact that the same set of templates was 
used to create the simulated data and to model them. In 
principle, we assumed that the grid of templates is repre- 
sentative of the real stellar population of the analysed spec- 
tra, but it is very possible that our library is incomplete. 
This is not the case in our simulations, by construction. To 
study the effects of this possible incomplete representation 
of the observational spectra by the template library, we cre- 
ated a new set of simulated data, corresponding to a flux 
level of 10 10~^® erg s"^ cm"^ A"\ that was then fltted 
using a much reduced version of the template library. This 
template comprises a grid of 12 SSPs, corresponding to 3 
ages (17 Gyrs, 900 Myrs and 25 Myrs), and 4 metallicities 
{Z =0.0001,0.004,0.02,0.05). The resuhs from this simula- 
tion are listed in Table A3, and shown in Figure A3. It shows 
that in the case of an incomplete coverage of the spectro- 
scopic parameters by the model template, these parameters 
are still well reconstructed, albeit with a lower accuracy (as 
expected). 

Finally, we tested how well the code is able to reproduce 
a pure SSP, instead of a mixture of them. With this test we 
try to overpass the intrinsic problems in the interpretation of 
the luminosity-weighted ages /metallicities in galaxies, and 
to test if we can recover the better understood properties of 
a pure single-burst population. The results from this simula- 
tion are listed in Table A3, and shown in Figure A4. Despite 
the discrete coverage of the input parameters, it is clear that 
the parameters are recovered just slightly worse than in the 
case of a multiple component stellar population. 

On the other hand, the relative difference between the 
original spectra and the recovered model is lower than ~ 15% 
for all the simulations, even in those ones for which the ac- 
curacy of the recovered parameters is deficient. These dif- 
ferences are of just 2-3% for the spectra with the highest 
signal-to- noise (>100). This basically means that the gen- 
eral shape of the spectra is better recovered than the param- 
eters that defines them. This is easy to understand, since 
these parameters are characterised by particular features in 
the spectra, rather than by their general shapes, which are 
more degenerated. 

Figure A5 shows the radial distribution of the average 
signal-to-noise ratio at ~ 5000 A of the spectra within the 
flux-calibrated mosaic of NGC 628, as described in Sec. 3. 
On the basis of this radial proflle and the results of the sim- 
ulations discussed before it is possible to recover the con- 
sidered parameters of the composite stellar populations for 
individual spectra only at radii less than ^ 10 arcsec, for the 
central kpc of the galaxy. The results are robust only for the 
integrated spectrum of the galaxy (with a S/N~400), and 
for radii lower than ~20 arcsec (with S/N~ 100). However, 
for the azimuthally averaged spectra, where the signal-to- 
noise increases by the root-mean square of the number of 
individual averaged spectra, the simulations indicate that 
the recovery of the considered parameters is robust even at 
radii as large as ~ 120 arcsec, above which both the S/N and 
the number of spectra with accurate spectrophotometry also 
drops (see Figure 4). The combination of both effects im- 
pose a limitation to the accuracy of the derived parameters 
at higher radii. 



